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Abstract 



This dissertation focuses on phenomenological studies for possible signals for super- 
symmetric events at the Large Hadron Collider (LHC). We have divided our en- 
deavours into three separate projects. First, considering that the branching fraction 
for the decays of gluinos to third generation squarks is expected to be enhanced in 
classes of supcrsymmetric models where either third generation fermions are lighter 
than other squarks, or models of mixed higgsino dark matter which are constructed 
in agreement with the measured density of cold dark matter (CDM), the gluino 
production in such scenarios at the LHC should be rich in top and bottom quark 
jets. Requiring 5-jets in addition to missing energy EJ^^^^ should, therefore, enhance 
the supersymmetry signal relative to Standard Model backgrounds. We quantify 
the increase in the supersymmetry reach of the LHC from 6-tagging in a variety of 
well-motivated models of supersymmetry. We also explore top-tagging at the LHC. 
Second, we explore the prospects for detecting the direct production of third gen- 
eration squarks in models with an inverted squark mass hierarchy. This is signalled 
by 6-jets -I- £'™'^^ events harder than in the Standard Model, but softer than those 
from the production of gluinos and heavier squarks. We find that these events can 
be readily separated from SM background (for third generation squark masses in the 
200 — 400 GeV range) , and the contamination from the much heavier gluinos and 
squarks although formidable can effectively be suppressed. Third, we attempt to 
extract model-independent information about neutralino properties from LHC data, 
assuming only the particle content of the MSSM and that all two-body neutralino 



decays are kinematically suppressed, with the neutrahno inclusive production yield- 
ing a sufficient cross section. We show that the Lorentz invariant dilepton mass 
distribution encodes clear information about the relative sign of the mass eigenval- 
ues of the parent and daughter neutralinos. We show that we can extract most 
neutralino mass matrix parameters if there is a double mass edge. 
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Chapter 1 

Introduction to Supersymmetry 



1.1 Introduction 

Supersymmetry (SUSY) [U [21 O HI [3] is one of the more peculiar theoretical dis- 
coveries in the history of physics, since, despite the enormous effort invested in its 
study (its discovery dates over 30 years ago), there is no experimental evidence of 
SUSY. 

Additionally, in the past decades, the Standard Model (SM) has been verified with 
great precision by numerous experiments. When discrepancies have been encoun- 
tered, these vanish with the increased precision in measurements, and the greater 
the precision that is achieved, the more precisely is the SM confirmed. From the 
viewpoint of precision measurements, there is little need for new physics beyond the 

sa£1 

The experimental evidence of the need for new physics beyond the SM comes from 
neutrino physics and observations supporting the existence of Dark Matter in the 
Universe without viable candidates in the SM. Also gravity interactions are not part 
of the SM. From the theoretical perspective there are also good reasons for going 
-"^Exceptions to this are the g — 2 tZj,^ experiment and the proton size anomaly [S] 
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beyond the SM: the solution of the hierarchy problem, or the desire to find a new 
unified theory, or a much simpler one, that offers an explanation of the symmetries, 
the spectrum, or the parameters of the SM. SUSY is one of the best candidates we 
have to this date of new physics beyond the SM. Besides offering a natural solution 
to the hierarchy problem, it allows for the unification of the gauge couplings. 
We must keep in mind that the hierarchy problem was not the primary motive be- 
hind the invention of SUSY in the 1970's. This is why it is surprising that although 
the initial SUSY models were quite different from the current Minimal Supersym- 
metric Standard Model (MSSM), with time this has become the principal candidate 
to succeed the SM. With the advent of the new generation of colliders, such as the 
Large Hadron Collider (LHC) coming online, a more definite test of SUSY is viable. 



1.2 SUSY Theory 

The construction of a SUSY theory would have as its underlying algebraic structure 
that of a graded Lie Algebra (gLA). These are extensions of the Lie Algebras, in 
which a distinction is established between elements of odd and even nature. Those of 
even nature obey commutation rules (Lie Algebra), while those of odd nature obey 
anti-commutation rules amongst them, and commutation rules with the even ones, 
i.e. the elements of odd nature constitute a representation of the gLA, so that for 
Ajn and Qa being the even and odd elements of this gLA respectively, we would have 



[-^m,Qa]= ^maQlS 
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{Qa,Ql^}= F^pA^ (1.1) 

where the repeated indices on the right-hand side are summed over. In the context 
of the extensions of the Poincare group by one self-conjugate spinor charge Q the 
even generators are the generators of the Poincare group and the odd generator is 
the generator of SUSY. A supersymmetry transformation will turn a bosonic state 
into a fermionic state, and viceversa, with the generator given by the operator Q as 
an anticommuting spinor, so that, 

fermion -f-^ boson 

An extension of the Coleman-Mandula theorem |10j . by Haag-Lopuszanski-Sohnius 
[TT] . restricts the possible supersymmetries acceptable in a Quantum Field Theory 
(QFT) with interactions. Only theories with one spinorial charge Qa, known as = 
1 SUSY, allow for chiral fermions, i.e., fermions whose left-handed and right-handed 
pieces transform differently under symmetry transformations, theories crucial for 
phenomenology. For this reason, we restrict our focus to A^ = 1 SUSY. We can then 
write the algebra as, 

{Qo.,Qp]= 2(7^)a/3PM (1-2) 

The algebra closes to yield the generators of the Poincare group, so these 

show that supersymmetry is a spacetime symmetry. The irreducible representations 

of the SUSY Algebra are labeled supermultiplets, each containing both bosonic 

and fermionic states having the exact same number of degrees of freedom. The 

supermultiplets of the SUSY Algebra utilized in the construction of the MSSM are: 
^the M^^'^ are the generators of Lorentz transformations 
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• ( <l>, ^' ) chiral superfield (or scalar , or matter) consisting of one Weyl fermion 
(uf = 2) and two real scalars (n^ = 2x1). 

• ( V^,\ ) vector superfield (or gauge) consisting of one spin-1 massless boson 
(rift = 2) and one Weyl fermion (nj = 2). 

The operator commutes with all generators, so that all particles occurring in a 
supermultiplet will have the same eigenvalues of P^, and therefore the same mass. 
The supersymmetry generators also commute with the genrators of gauge transfor- 
mations, which means that members of a supermultiplet belong to the same repre- 
sentation of the gauge group, thus having the same electric charge, weak isospin, 
and color degrees of freedom. 

The superfield formalism [12] provides a convenient tool for studying supersymmet- 
ric theories including studying the multiplet structure of these theories, the unitary 
supersymmetric representations of the particle states, the construction of super- 
symmetric invariants, amongst others. Most importantly, it provides a recipe for 
the construction of a Lagrangian density of a supersymmetric Yang-Mills theory 
in terms of the ordinary boson and fermion fields of QFT. Assume that the chiral 
supermultiplets (cpi^ipi) transform under a gauge group representations and that 
{V^,Xa) are the gauge supermultiplets, with a as the gauge group index. Then the 
Lagrangian density can be written as. 



= + + Cy + Cs + Cmg + ^XG + ^GG 



(1.3) 



where Ck contains the kinetic term. 



\v^,Vr + \\''^\a + h.c. (1.4) 
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Cm\ gives the interactions of gauginos with the scalars and fermions of the chiral 
multiplets 

Cmx = -ig V2{(P* {ta)T iPl, + h.c. (1.5) 
Cy yields the fermion mass terms and the Yukawa- type interactions, 

Cs includes interactions between scalar fields , known as F-terms and D-terms , 
which have an important role in the breaking of SUSY 

Cmg gives us the interactions between each particle and the gauge fields, 

Cmg = -g r {ta)i, - ig 4>*i {ta^ + 

+5V;y^^ 0* {tah)ij (Pj + h.c. (1.8) 

CxG includes the interactions between the gauge fields and the gauginos, 

CxG = igfabcX''r>'v^ + h.c. (1.9) 

and Cgg contains the self-interactions of the gauge fields where the are the 
matrices of the Lie Algebra associated with the gauge group, 

JC-GG = -g fL d^V: - f^, fade V^V^F^'^y- (1.10) 

which satisfy 

[ta,tb]=irabtc (1.11) 

The complete lagrangian may be written in a more compact form if covariant deriva- 
tives are used, as illustrated in eqs. (6.44) and (6.45a-d) in [1]. 



The model is completely specified once the superpotential is specified. For 

chiral superfields renormalizability requires W to be the most general gauge invari- 
ant function that is at most cubic in the fields, i.e. 

= 5Z ^i^i ^ ^ "^^-^ ^^^^ Kjkhh^k (1-12) 

1.3 Spectrum of the MSSM 

To create a SUSY version of the SM we must pair up the SM fields in supermultiplets 
and introduce the SUSY partners of all the SM fields. For the 12 gauge bosons of 
SU{3)c X SU{2)l X U{1)y there are no available partners so we must introduce 
12 fermions: 8 gluinos {g), 1 bino (Aq), and 3 winos iWj), whose definition is a 
generalization of the one in the SM. For the chiral fermions of the SM we need 
new complex scalar partners known as squarks and sleptons. One Higgs doublet 
is not enough to provide mass to both of the charge (-1/3 , 2/3 ) quarks and 
satisfy invariance under SUSY transformations, so we must introduce an additional 
Higgs doublet, defined in Table |1.1[ This is also, what is required to cancel the 
chiral anomaly that would otherwise arise. The resulting spectrum, as the simplest 
SUSY generalization of the Glashow-Weinberg-Salam model then consists of the 
fields listed in Table 11.11 

1.4 A SUSY Toy Model 

As an illustration consider a field theory [T, 13] with Lagrangian given by, 

— C'kin ~\~ C'mass 
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Field 


SU{3)c X SUi2)L X U{1)y 




(1,2,-1) 




(1,1,2) 




(3,2,^) 




(3M,-i) 




(3M,i) 




(1,2,1) 




(1,2*,-1) 



Table 1.1: MSSM particle content. Only the first generation of matter particles is 
shown; the second and third generations are replicas of this. 
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with 

Ck^n = lid.Af + ^{d,Bf + liF' + G') (1.13) 

Cmass = -m[hp^l; -GA- FB] (1.14) 

where A and B are real scalar fields with mass dimension [A\ = [B] = 1, ^/j is a 4- 
component Major ana spinor field with mass dimension [-0] = 3/2, and F and G are 
also real scalar fields with dimension [F] = [G] = 2. Both F and G have no kinetic 
terms, so their equations of motion are algebraic and can be used to eliminate these 
fields from the Lagrangian. This yields, 

C = \{d^Af + \{d^Bf M- \^\A^ + B^) - \m'-^i^ (1.15) 

As a note, we see that the number of bosonic and fermionic degrees of freedom 
in the Lagrangian exactly balance: without the equations of motion, the four real 
components for the Majorana spinor field are balanced by the four real scalar fields. 
This is the Lagrangian for free fields A, B and tjj. After applying the respective 
equations of motion, their quanta correspond to two spin zero particles A and B 
and a self-conjugate, spin ^ particle, all with the same mass. We can add interactions 



to our Lagrangian (1.15), assuming renormalizability, 

Ant = -^^V^^ + ^^V^TbV' - gmV2AB^ - ^A{A^ - B^) 

-g^A^B'-^{A'-B^f (1.16) 



Note that the Lagrangian (1.16) has just one mass and one coupling parameter for 
all fields. 

It is important to note that in the SM radiative corrections to the mass of the 
fundamental scalars of the SM will be quadratically divergent. If the SM is coupled 
to new physics at a scale M, the quadratic divergence manifests itself as corrections 



that grow as M^, and so destabilize the weak scale if M is much larger than the Fermi 
scale, which is the natural scale of the SM. This is known as the hierarchy problem. 
The expectation is, that the new physics will have inherent a symmetry which will 
induce the cancellation of the quadratically divergent contributions. SUSY is just 
such a symmetry, and as long as the new energy scale is < 1 — 2 TeV, the fine-tuning 
is ameliorated. 

To illustrate the cancellation of quadratic divergences, we can use our toy theory 
and show the cancellation of the quadratically divergent contributions that would 
destabilize the scalar sector. Before doing so, we need to evaluate the following 
quadratically divergent integral up to some momentum cutoff value A, 

^id = [ TTT-fj -2 2 I ■ (1-17) 

J {27rp — + te 

Noting the poles of the integrand and choosing an appropiate contour, and setting 
limits of integration, 

dS i 1 



/ 



{2tiY - m? + ie in'^ Jo y^q'^ + m'^ 

lad ~ -Kl^^ - rn^lnl — ) + const x m^l (1-18) 



Proceeding as in [T] , we consider first the one-point function of the field A to first 



order in the coupling g. The relevant interaction Hamiltonian from (1.16) is 



Hint = -Hint 3 ^^V'V' + ^"^^^^ + ^""^^ (^-^^^ 

If we expand the matrix element (r2|T^(x)|r2), where |r2) is the ground state of the 
interacting theory, perturbatively to order G, we get. 



V2 



1 dSB^{x - y)\{-\)TrSF{y - y) + mi^f (y - y) + ?>mD^{y - y)\ (1.20) 
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where Dp is the Fourier transform of the scalar field propagator in momentum space 
given by, 



Df{x - y) 



-iq-(x~y)_ 



(2vr)4 

The factor in square brackets above is proportional to, 



^ — + ie 



(1.21) 



Tr 
4m^ 



d'^p 



d'^p 
p^ - m'^ 



m 



m 



d'^p 



p^ — 

d^p 
p^ — 



3m 



3m 



d^p 



p^ — m\ 

d'^p 
p^ — m^ 



(1.22) 



where m^,mA,i^B ai'e exactly the same as mass parameter m in the trilinear scalar 



couplings in eq. (1.16). Since these masses are exactly equal in a supersymmetric 



theory, the three contributions in (1.22) add to zero. So although each contribution 
is separately quadratically divergent, the divergence due to the fermionic term can- 
cels the sum of the divergences from the bosonic terms. In order for this to happen 



it is necessary that the couplings are exactly those in (1.16). Also, the quadratic 



divergence in (1.22) is independent of mA and m^, however the fermion mass must 
be equal to m. 

If we look at the lowest order quadratic divergences in the two-point function of 
{Q\TA{x)A{y)\^), once again the quadratic divergences cancel out between 
fermionic and bosonic contributions, and the cancellation occurs for all values pf 



particle masses. It is again crucial that couplings are as in (1.16). We thus see 
that as long as the dimensionless couplings are as given by supersymmetry, the 
quadratic divergences cancel even if supersymmetry is broken by scalar masses dif- 
ferent from fermion masses. This is an example of Soft Supersymmetry Breaking 
(SSB) discussed in the next section. 
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1.5 SUSY Breaking 

If SUSY were an exact symmetry of nature, SUSY particles would have the same 
mass as their SM partners. This is not so, otherwise discovery of the spartners of 
the known particles should have been possible at the accelerators available prior 
to the LHC operation. Therefore, SUSY is a broken symmetry at the Fermi scale. 
Having a simple model to describe the breaking of SUSY which we could connect to 
the MSSM would make our endeavours much easier, but unfortunately, such models 
are far from being simple. As previously mentioned during the discussion of our toy 
model, we can add to JCsusy terms which violate supersymmetry but which are of 
little importance at high energies. The complete list of possible terms, all of which 
are required to have mass dimension < 4, which may be added to CguSY without 
altering the cancellation of quadratic divergences in the radiative correction to the 
SM Higgs mass is as follows: 



• mass scalar terms: 4't^3 ' ^i^j 

• trilinear scalar interactions: aijk(pi(pj4'kj Cijk4'*i4'o'Pk and their h.c. 

• gaugino masses: \Mi\i\i + h.c. 

• linear terms: Cj(^i 

These all are known as Soft-SUSY-Breaking terms (SSB). We then distinguish two 
separate components in the complete Lagrangian density: 

^ = ^SUSY + ^SOFT 
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CsoFT terms parameterize the fundamental mechanisms of SUSY breaking and in- 
clude the majority of the parameters appearing in the Lagrangian. This complete 
Lagrangian is what we have introduced before as the MSSM. Any sensible phe- 
nomenological study is impaired by the very huge number of parameters which are 
present in CsoFT- It would be helpful to have a theory capable of predicting the 
soft parameters, and effectively there are several such theories. 
Two fundamental mechanisms exist which attempt to explain how MSSM super- 
partners acquire their masses. Common to both of them is the existence of a hidden 
sector responsible for SUSY breaking and an interaction responsible for transmitting 
the breaking to the visible sector where we find the particles which constitute the 
MSSM. These fundamental models are known as, 

1. Gauge Mediated SUSY Breaking (GMSB): where the transmission interaction 
is the same gauge interaction of the SM jl5j . 

2. Gravity Mediated SUSY Breaking (SUGRA): where it is gravity which acts 
as the messenger for SUSY breaking |16j . 

In these models the parameters of the MSSM are determined in terms of a handful 
of parameters at specified high energy scales. A consequence of this is that the Higgs 
mass parameters acquire negative values and produce the ElectroWeak Symmetry 
Breaking (EWSB). We see then that EWSB is intimately related to SUSY breaking. 
The minimal SUGRA (mSUGRA) model has been extensively studied phenomeno- 
logically, as well as the GMSB model. Within mSUGRA, the soft parameters acquire 
a simple structure at the unification scale (or Planck scale) in which. 
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1. Scalar masses are universal (diagonal) 

•mq = m% = ml = ml = rn% = m% = ml 

2. Gaugino masses are universal 

Ml = M2 = M3 = mi 
2 

3. The cube terms of the soft potential are proportional to the superpotential 
Yukawas 

[a/] = Mh}] 

where ^0 is a common parameter. At low energies, parameters are determined 
by the renormalization group equations (RGE) from their high energy values. For 
mSUGRA, the model is fixed by the 18 parameters from the SM with five additional 
parameters 

(mo,mi/2, Ao,5o,/i/) 

Radiative EWSB (Electroweak Symmetry Breaking) determines /x^ and it is tradi- 
tional to eliminate Bq in favor of tan/3 leaving the often used parameter set 

(mo, mi/2, ^0, ton/3, sgn{ij)) 

If the sparticle masses are ~ 10^ — 10^ GeV, then extrapolating to high energies 
by way of the RGB's, with the three SM gauge couplings measured at the weak 
scale, these very nearly meet at a point under MSSM evolution, suggesting physics 
at scales Mqut ~ 2 x 10^^ GeV is described by a SUSY GUT. Below Mqut, the 
correct effective field theory is provided by the MSSM. 
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We would like to introduce a symmetry which acts differently on the component 
fields of the superfields so as to distinguish the SM particles from their superpart- 
ners. This is known as i?-symmetry, consisting of a U{1) phase applied to the 
components of the superfields. When the phase is restricted to the value vr, the 
i?-symmetry is referred to as i?-parity, with the phase being either +1 or —1. All 
SM particles are even under i?-parity while superpartners are odd under i?-parity. 
If we now restrict CsusY to be invariant under i?-parity, a consequence will be that 
there cannot be interactions coupling a single superpartner to two SM particles. 
This implies that all superpartners will ultimately decay to an sparticle, the lightest 
supersymmetric particle or LSP, which will be stable. Then the Universe must be 
filled with these sparticles, and from experiments on the charge-to-mass ratio of 
matter, the possibility of it being electrically charged has been ruled out. Thus, the 
LSP is electrically neutral. The LSP would be a viable candidate for Dark Matter. 
In our work, we assume i?-parity invariance. A formula to calculate i?-parity is 
given by 

R = (_l)3(B-i)+2s 

where B, L are the baryon (lepton) quantum numbers and s is the spin. The factor 
(— 1)^* guarantees that particles and their superpartners will have opposite i?-parity. 
Imposing i?-parity invariance eliminates all baryon (lepton) violating interactions, 
if interactions are renormalizable. An important phenomenological consequence of 
i?-parity invariance is that sparticles can be produced only in pairs at colliders, and 
must decay to SM particles plus an odd number of sparticles. 

Summarizing, the MSSM is a quantum field theory with supersymmetry relating 
bosonic and fermionic degrees of freedom. Providing us with a solution to the SM 
fine-tuning problem by eliminating the quadratic dependence on the cut-off scale A, 
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it is perturbatively calculable for all energies up to Mpia„ck or Mqut scales without 
requiring huge fine tuning. It is less UV divergent than corresponding non-SUSY 
theory due to cancellation of the leading quadratic divergence of fermionic loops 
with those of bosonic loops. When elevated to a local supersymmetry gravity is 
automatically introduced leading gravitational interactions along with strong, weak 
and electromagnetic interactions in an effective field theory. With conservation of 
i?-parity it includes a stable massive particle which is usually electrically and color 
neutral, providing us with an excellent candidate for the observed cold dark matter 
in the Universe. 
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Chapter 2 

SUSY Phenomenology at the 
LHC 

2.1 Large Hadron Collider 
2.1.1 Overview 

Located in Geneva, Switzerland, the Large Hadron Collider LHC is a pp collider, 
built to operate initially at a CM. energy of y/s = 7 TeV, with plans to eventually 
reach its design energy of ^/s = 14 TeV. At these energy scales, the LHC is capable 
of creating the conditions essential for discovery of new physics at the weak scale 

[nms]. 

The discovery of a single SUSY particle would be as groundbreaking as that of find- 
ing the elusive Higgs boson. The available energy would facilitate the discovery of a 
superpartner particle, allowing these hypothetical particles to manifest themselves 
over the background from SM sources. Such an event holds promise of explaining 
one of the most persistent mysteries in physics and astronomy, the existence of dark 
matter, first theorized in the 1930's. 

The LHC work may also reveal the existence of additional dimensions of space, if 
nature really has hidden dimensions, over the known 3 + 1, allowing for a struc- 
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ture of the Universe more complex than our current knowledge describes. The LHC 
has a length of 26.659km, with detectors located at several points along the cir- 
cumference. The proton beams are obtained by ionizing the Hydrogen gas travel 
in opposite directions while being accelerated to speeds very close to the speed of 
light. This happens in a successive series of smaller accelerators, before the beam 
is finally injected into the LHC for the last stage of acceleration. Here, powerful 
frequency devices provide a kick to the particles each time they pass by. There are 
over 10-thousand superconducting electromagnets, supercooled by liquid Helium to 
1.9K. 

The detectors are: 

• Compact Muon Solenoid (CMS) [HEI]. 

• A Toroidal LHC Apparatus (ATLAS) [201 ESI [23] . 

• A Large Ion Collider Experiment (ALICE) [Ml ES [26]. 

• LHCb[27l[28]. 

CMS and ATLAS are general purpose detectors and the analysis of data of these 
experiments will be of primary interest to us in this dissertation. 
With the LHC using approximately 2x 10^^ grams of Hydrogen per day, the ultimate 
collisions emerge after a succession of processes: 

1. Protons are given an initial boost in the small linear accelerator known as 
Linac2. to about 0.314c. 

2. Protons then move into CERN's old circular accelerator, the Proton Synchro- 
ton (PS)Booster, where they can boost their speed every lap until reaching 
speeds of about 0.916c. 
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3. The next boost happens at the Proton Synchroton, to about 0.9993c. 

4. Protons are now funneled into the Super Proton Synchroton (SPS), where in 
1983 both the W and the Z were first detected. Here protons reach ~ 0.99998c 
which is equivalent to a CM. energy of about 450GeV. 

5. Protons are led into the LHC, where for a CM. energy of 7 TeV at current 
operation, the speeds are ^ 0.999999991c . One beam consists of some 2808 
bunches, with w 10^ protons in each bunch. 

6. Collisions now occur at ATLAS (point 1), CMS (point 5), LHCb (point 8), 
and ALICE (point 2). 

The ability of a detector to find and measure particle momenta with high accuracy is 
propotional to the strength of the magnetic field Bx the distance travelled inside the 



detector. For the CMS design, as shown in Fig. 2.1 the choice was to build a compact 
instrument offering a relatively short path for the muons inside the detector, but 
using a high magnetic field B ~ ATesla. 



The ATLAS detector, as shown in Fig. 2.2 was designed with the alternate choice, 



a bigger instrument offering a larger path, but using a smaller magnetic field B ~ 
ITesla, thus achieving the same capability. The pre-assembled piece of the CMS 
detector, containing the giant electromagnet weighs in as much as 5 Boeing 747 
airplanes. 



2.1.2 Techniques for LHC Searches 

From the billions of pp collisions, some 10-15 petabytes (lO^^bytes) of data are gen- 
erated per year. In the CMS, from this huge amount of collisions, maybe 1 in 10^ is 
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Figure 2.1: CMS detector at the LHC 

of interest, 300 are permanently recorded for complete reconstruction and analysis, 
and 1 is placed on screen every second. 

Four kinds of emissions follow proton collisions in the LHC: 

1. Jets: streams of quarks and gluons that ultimately convert to hadrons emanat- 
ing at various angles, depending on the energy and type of reactions produced 
in the collision. 

2. Discrete emissions of isolated leptons : e, //, r. 

3. Missing transverse energy E™*^*^: energy of particles that are undetected and 
moving in directions that have a transverse component to the direction of the 
colliding beam of protons. 

4. Photon emission. 
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Figure 2.2: ATLAS detector at the LHC 

The total transverse momentum of the final products in the center of mass frame of 
the collision should be zero. Then the difference between the measured amount and 
zero yields the i?™'^^(or alternatively, -^rp). In hadronic collisions, the partons which 
participate in the hard process carry a fraction of the beam energy. The remnants 
of the beam associated with the remaining partons mostly escape undetected in the 
beam pipe. Thus, only conservation of the momentum in the direction transverse 
to the beams is relevant, making the missing transverse energy i?™'^^ the important 
quantity, rather than the total missing energy. 

February 2010 provided the first report of a collision taking place inside the 
CMS detector at the LHC late in 2009. On March 19, 2010 LHC reached its target 
energy for the next two years : 3.5 TeV per beam, 7 TeV total. After this period 
the LHC will undergo a year of maintainance, following which it expects to operate 
at its design energy of 14 TeV for both beams. After more than a year of operation, 
its total integrated luminosity for 2010-11 is 3.3 fb~^, as of this writing. Recently 
recorded events are viewable at [29]. 
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The ATLAS and CMS have searched for an excess of events above SM expectations 
in channels that would be populated by the production of gluinos and squarks at 
the LHC. Unfortunately no such excess was found. This has then been translated to 
upper limits on the cross section for particle production. This exclusion is illustrated 
within the mSUGRA model, introduced in Chapt.[T} in Fig 2.3 for ATLAS data and 



in Fig 2.3 for CMS data where a composite is made for an integrated luminosity of 
about 1 fb^^ that was analyzed in Summer, 2011. 
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Figure 2.3: Exclusion region in the mSUGRA/CMSSM {mQ,mi/2) plane for tan (3 
10, ^0 = and > for an integrated luminosity of 1.04 fb~^ at ATLAS [30] 



2.2 SUSY Event Simulation 

The key link between theoretical predictions of SUSY or other new physics and the 
actual experimental observations of particle tracks and calorimeter depositions in 
collider detectors is provided through event generator programs. These allow us 
to compute how a theory would manifest in actual collider experiments. With the 
LHC operating at energy scales of ^/s = 7—14 TeV, this should be enough to pro- 
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Figure 2.4: Exclusion region in the mSUGRA/CMSSM (mo, mi/2) plane for tan /3 = 
10, Aq = and fi > for an integrated luminosity of 1.04 fb~^ at CMS. 



duce superpartners and provide evidence for viable particle models, such as weak 
scale SUSY. It is possible that discovery of new physics beyond the SM can result 
from indirect searches, but for SUSY at the weak scale it is widely accepted that 
evidence for it will come from direct creation of supersymmetric matter in colliding 
beam experiments., and the detailed analysis of the resultant scattering events. [1] 



2.2.1 Event generator 

Different models describing supersymmetry exist, which are used to predict sparticle 
production rates and subsequent decay patterns into final states of SM fundamental 
particles. Some of these, such as quarks and gluons, cannot be detected directly in 
a collider detector. The detectors will measure tracks and momenta of quasi-stable 
charged particles bending in a magnetic field, in addition to energy deposited in 
calorimeter cells by charged leptons, hadrons, and photons. We realize then that a 
gap exists between the predictions of SUSY theories for final states of fundamental 
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particles, and the actual detection at the experimental level. It is the existence of 
this gap that necessitates the development of event generator programs[31j. Cur- 
rently available general purpose event generator programs that incorporate SUSY 
include IS A JET [31], PYTHIA 33J, HERWIG [IJ, SUSYGen [33], and SHERPA 
Once a SUSY theory and collider type are specified, the event generator program 
produces a full simulation of types of scattering events to be expected. The final 
states are completely specified including detailed kinematics of each particle in the 
event. 

The present work studies prospects for physics at the LHC operating at its design 
energy of 14 TeV in the center of mass. Then, for an input of an MSSM set of 
parameters, the generator, in our case ISAJET version 7.74, generates the sparticle 
pair production events according to the ratio of their production cross sections. The 
sparticles will then undergo a decay into a partonic final state, according to branch- 
ing ratios specified by the model. Then the partonic state final state is converted to 
one composed of particles which are detected experimentally. 

The fraction of the hadron's longitudinal momentum carried by the initial hard 
scattering partons is unknowable, so there is an irreducible uncertainty in the longi- 
tudinal boost of the center-of-momentum frame for the colliding system. However, 
by forming the vector sum of all the energy deposited in the transverse direction, 
we get an important quantity called missing transverse energy, ^rp, as mentioned 
previously. A certain amount of-^rp is due to jet and lepton mismeasurement from 
imperfect energy resolution, particles going into cracks in the detector, and other 
non-physics causes. A large .^rp, however, generally indicates the production of one 
or more high-energy weakly-interacting particles that escape the experimental appa- 
ratus without depositing energy. In the Standard Model these would be neutrinos. 
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In SUSY searches, large -^rp is the signature of escaping LSPs. Indeed, since a pair 
of LSPs are always produced in a SUSY reaction where ii-parity is conserved, large 
is the hallmark of a supersymmetric reaction. 



2.2.2 Detector simulation 

As each event is generated, it is processed through a toy detector simulator consisting 
of the following elements (where rj is the pseudorapidity, (j) is the azimuthal angle, 



and AR = ^{Arf)"^ + (A^)^). The toy detector captures the salient features of the 
LHC detectors. 

• Calorimeter simulator: We implement a toy calorimeter based on the ISAJET 
CALSIM subroutine. The segmentation is Ar/ x A(/> = 0.05 x 0.05 extending 
to a rapidity of \ri\ = 5. There is a hadronic calorimeter, into which hadrons 



deposit their energy with a resolution given in Table 2.1 for different ranges of 



1 7/ 1, and an electromagnetic calorimeter which captures electrons and photons 



with resolution also listed in Table 2.1 We do not attempt to simulate effects 



of cracks or dead regions that are specific to particular detectors. 

Isolated lepton identification: We sum the hadronic transverse energy in a 
cone of Ai? < 0.3 around each lepton. If this hadronic energy is less than 50% 
of the lepton's transverse energy, then the lepton is declared isolated. The p^^ 



thresholds for isolated leptons are given in Table 2.1 for each case 



Jet identification: We use ISAJET's GETJET jet-finding algorithm. Jets are 
defined as hadronic clusters with total Er^ > 50 GeV falling within a cone of 
radius AR < 0.5 and subject to |?7| < 3. We do not correct jet energy. 
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Silicon vertex detector (SVX): We simulate a SVX detector for tagging 6-jets. 
We identify each weakly-decaying B hadron in an event with Ej, > 15 GeV 
and 1 7/1 < 3. If Ai?(i?, jet) < 0.5 for some jet then that jet is tagged as a 6-jet. 
with an efficiency of 50% [37] at the LHC design luminosity of 100 fb^-^/y, and 
assume that gluon and light quark jets can be rejected as b jets by a factor 
Rb = 150 (50) if Et < 100 GeV {Et > 250 GeV) and a linear interpolation in 
between [37] . 



In Table |2.1| we summarize the basic parameters used to define jets, b-jets, and 
isolated leptons. 



calorimeter \r]\ 
cell size 

hadronic resolution 

\v\ < 3 
3< \ri\ < 5 

e 

EM resolution 

Jets are 

Et >50 GeV 

|7?| <3.0 

B-Jets are 

Et > 50 GeV 

B-hadron 

Isolated Leptons 

in LHC reach studies 

in LHC dilepton studies 



< 5 

Ar]xA(j) = 0.05 x 0.05 



O.5/V^0O.O3 
1/V^0O.O7 
addition in quadrature 
O.l/^/^0O.Ol 

hadronic clusters 
AR = VA??2 + A(/)2 = 0.5 

tagged at 50% eff 
\r]\ < 3 
PT > 15 GeV 

PT > 10 GeV 
PT> Q GeV 



Table 2.1: simulation initial set of parameters defining jets bjets and isolated leptons 
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2.3 Phenomenological goals 

Given the opportunity to make a wish hst of measurements that may shed hght on 
the properties of the sparticles, which would estabhsh a roadmap to the high energy 
scale physics responsible for the breaking of SUSY, one would set apart as being of 
upmost importance the following: 

• The discovery of sparticles, whose properties would be indicative of the validity 
of existing SUSY models. 

• The measurement of the masses of as many superpartners as possible. 

• The relevance of signals from third generation squarks, possibly shedding some 
light on an Inverted Mass Hierarchy (IMH) in SUSY, where the sparticles 
belonging to the third generation have masses lower than the corresponding 
ones in the first/second generation. 

• the determination of the parameters involved in higgsino-gaugino mixing 

In view of the above as motivation for this dissertation, we have divided our endeav- 
ours into three separate projects. 

1. Heavy-flavor tagging and the SUSY reach at the LHC: 

Considering that the branching fraction for the decays of gluinos to third 
generation squarks is expected to be enhanced in classes of SUSY models 
where either third generation squarks are lighter than other squarks, or models 
of mixed higgsino dark matter which are constructed in agreement with the 
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measured density of cold dark matter (CDM), the gluino production in such 
scenarios at the LHC should be rich in top and bottom quark jets. Requiring 
6-jets [38] in addition to Ej^^^^ should, therefore, enhance the supersymmetry 
signal relative to Standard Model backgrounds from V + jet, VV and QCD 
backgrounds {V = W,Z). We quantify the increase in the supersymmetry 
reach of the LHC from 6-tagging in a variety of well-motivated models of 
supersymmetry. We also explore "top-tagging" at the LHC. We find that 
while the efficiency for this turns out to be too low to give an increase in 
reach beyond that obtained via 6-tagging, top-tagging can indeed provide a 
confirmatory signal if gluinos are not too heavy [39] . 

2. Signals for light third generation squarks (stops) at the LHC: 

We explore the prospects for detecting the direct production of third gen- 
eration squarks in models with an inverted squark mass hierarchy. This is 
signalled by 6-jets + Ejl^^^ events harder than in the Standard Model, but 
softer than those from the production of gluinos and heavier squarks. We find 
that these events can be readily separated from SM background (for third 
generation squark masses ^ 500 GeV), and the contamination from the much 
heavier gluinos and squarks although formidable can effectively be suppressed 

m- 

3. Neutralino mass reconstruction and MSSM parameter deter- 
mination: 

We attempt to extract model-independent information about neutralino prop- 
erties from LHC data, assuming only the particle content of the MSSM and 
that all two-body neutralino decays are kinematically suppressed, with the 
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neutralino inclusive production yielding a sufficient cross section. We show 
that the Lorentz invariant dilepton mass distribution encodes clear informa- 
tion about the relative sign of the mass eigenvalues of the parent and daughter 
neutralinos. We attempt to answer question as to whether from the dilepton 
distribution we can establish if the decay is the result of a virtual Z-boson or 
a virtual slepton Il^h exchange. We attempt to extract information as to the 
values of the MSSM parameters that determine the mass of the neutralinos. 

The first two items listed above make up the contents of Chap.jS} while the third item 
is elaborated upon in Chap.|4j We end in Chap. [5] with a brief outlook for the future. 
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Chapter 3 



Heavy Flavor Tagging and the 
LHC Reach 



3.1 Introduction 



In the previous chapter, in Sec. 2.3 we presented three projects which comprise the 



focus of the work presented in this dissertation. Now, in this chapter, we detail our 



work with respect to the first two of these projects mentioned in Sec. 2.3 whose 
primary goal is to provide answers to the following two questions: 



• Ql: What information obtained from LHC events can we use to develop tech- 
niques which will allow us to extend the SUSY reach projections using Ej^^^^^ 
jets and leptons, at the LHC? Extending the SUSY reach encompasses ex- 
tending the region of MSSM parameter space where SUSY signals can be 
distinguished from SM background events, based on specific observability re- 
quirements. 

Additionally, we would like these techniques not to be restricted to a specific 
SUSY model, but rather to be applicable to classes of models. We consider 
models that: first, they conform to Dark Matter constraints, specifically CDM, 
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which we describe below. Second, we study models that exhibit an Inverted 
Mass Hierarchy (IMH), whereby the third generation scalar sparticles are 
lighter than their corresponding 1** and 2"*^ generation counterparts. IMH 
models are interesting because serves third generation sfermion mass parame- 
ters are driven to sub-TeV values, leaving first and second generation scalars as 
heavy as 2-3 TeV. The multi-TeV values of first and second generation scalar 
masses ameliorate the SUSY CP and fiavour problems without destroying the 
SUSY resolution of the gauge hierarchy problem, since the fields with sub- 
stantial direct couplings to the Higgs sector (gauginos and third generation 
scalars) have masses below the TeV scale. 

• Q2: Can we use these same techniques to isolate signals corresponding to direct 
production of third generation sparticles, not only from SM background events, 
but Additionally from other SUSY events which would now be considered an 
added contamination to the background? This would provide unequivocal 
evidence for the production of third generation squarks. 

It has been shown that both squark and gluino masses are smaller than 2 — 3 TeV 
and their production will be observable above SM backgrounds via signals consisting 
of multi-jet plus multi-lepton events with large amounts of Ej^^^^ carried off by the 
escaping LSPs. We will assume that the lightest neutralino is the LSP as is the case 
in many models. Remarkably, SUSY models with a stable neutralino LSP naturally 
lead to the right magnitude for the measured relic density of thermally produced 
cold dark matter [H], if superpartner masses are ~ 100 GeV. 

Defining 0, as the total matter/energy density of the Universe as a fraction of the 
critical closure density pc ~ 1.88 x 10^'^^ h^gcm~^ where h is the Hubble parameter 
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in units of WQkm/ sec/Mpc, then the component arising from non-relativistic and 
non-radiating matter is labeled as ilcDM with an inferred value of |43j . 

Qcdm/i^ = 0.1120 ±0.0056 , {2a) (3.1) 



Assuming thermal production and standard Big Bang cosmology, the upper limit 



from (3.1) provides a stringent constraint on any theory with stable weakly inter- 
acting particles, in particular on weak scale SUSY theories. Since the dark matter 
may well consist of several components, the contribution from any single component 
may well not saturate the observed value, so that strictly speaking the relic density 
measurement serves as an upper bound, 

< 0.12 , (3.2) 

on the relic density of neutralinos, or for that matter, on the density of any other 
stable particle. 



Direct searches for charged sparticles at LEP 2 have resulted in lower limits 
of about 100 GeV on chargino and selectron masses, and slightly lower on the masses 
of smuons and staus |45) . Since neutralinos can annihilate via t-channel sfermion ex- 
change, the measured value of the relic density, on the other hand, favours sfermions 
lighter than about 100 GeV, resulting in some tension with the LEP 2 bounds. In 
many constrained models where all sparticle masses and couplings are fixed by just 
a few parameters, such light sparticles often also lead to measurable deviations in 
other observables, and hence are disfavoured. If the SUSY mass scale is raised 
to avoid these constraints, the annihilation cross-section which is proportional to 
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-rJ^ — is correspondingly reduced, and the neutralino relic density turns out to be 
too large. One way to fix this is by invoking non-thermal relics or non-standard 
cosmology to dilute the relic density. However, it seems much more economical to 
invoke SUSY mechanisms that enhance the neutralino annihilation rate to bring 



their thermal relic density in line with (3.2) 



The primary reason for the low neutralino annihilation rate lies in the fact 
that the LSP is dominantly a bino in many models with assumed gaugino mass 
unification, where the bino and wino masses are related by Mi ~ The annihi- 

lation of bino pairs to gauge bosons is forbidden because SU (2) x U{1) precludes the 
couplings of binos to the gaugino-gauge boson system, while annihilation to fermions 
may be suppressed by large sfermion masses and the relatively small hypercharge 
coupling. Finally, annihilation to Higgs boson pairs is suppressed by the (usually 
large) higgsino mass, as well as by the small hypercharge gauge coupling. This then 
suggests several ways in which the neutralino annihilation rate may be enhanced to 



bring their thermal relic density in accord with (3.2). 



We can arrange the mass of a charged or coloured sparticle to be close to that 
of the LSP. Since these coloured/charged sparticles can annihilate efficiently, 
interactions between them and the neutralino which maintain thermal equilib- 
rium will necessarily also reduce the neutralino relic density [16]. Within the 
mSUGRA model, the co-annihilating sparticle is usually either the scalar tau 
|47j or the scalar top [48j, but different choices are possible in other models. 

We can arrange 2m ;y ~ ~ mn, so that neutralino annihilation is res- 
onantly enhanced through s-channel heavy Higgs boson exchange |49| . The 
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large widths of A and H together with the thermal motion of the LSPs in 
the early universe then enhances the annihilation cross section over a con- 
siderable range of parameters. Within the mSUGRA model, this is possible 
only if tan f3 is very large. However, in models with non-universal Higgs mass 
(NUHM) parameters, where the Higgs scalar mass parameters do not unify 



with matter scalar parameters as in mSUGRA \50\ I51|. agreement with (3.2) 
may be obtained via resonant A/H annihilation for any value of tan/3. We 
mention that resonantly enhanced annihilation may also occur via h exchange, 
albeit for a much smaller range of parameters [52] . 

It is also possible to obtain an enhanced neutralino annihilation rate if the 
light top squark, ti, is relatively light so that neutralinos efficiently annihilate 
via ZiZi —7- tt [53], or in NUHM models via ZiZi — )■ uu or cc, via t-channel 
top- or right-squark exchanges, respectively |51j . 



Instead of adjusting sparticle masses, we can also adjust the composition of 
the neutralino. More specifically, 

• We can increase the higgsino content of the neutralino so that its couplings to 
the gaugino-gauge boson pairs are increased, leading to mixed higgsino dark 
matter (MHDM). Within the mSUGRA framework, we can only do so in the 
so-called hyperbolic branch/focus point (HB/FP) region where mo takes on 
multi-TeV values |54j , but in NUHM models this is possible for aZZ values of uiq 
|51j . The higgsino content may also be increased by relaxing the assumed high 
scale universality between gaugino masses. The usually assumed universality 
of gaugino masses follows if the auxiliary field that breaks supersymmetry 
does not break the underlying grand unification symmetry; if this is not the 
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case, non- universal gaugino masses can result. It has been shown that if the 
GUT scale gluino mass is smaller than the other gaugino masses, fnjj^ does 
not run as negative as usual, yielding a smaller value of fj?, resulting in an 
increased higgsino content of Zi [55j. This has been dubbed as low M3 dark 
matter (LM3DM) . It has been pointed out |56] that increasing the GUT scale 
wino mass parameter from its unified value also results in a low value of |^|, 



resulting in consistency with (3.2) via MHDM. 



• Finally, depending on the gauge transformation property of the SUSY breaking 
auxiliary field, it may also be possible to enhance the wino content of the 
neutralino leading to mixed wino dark matter (MWDM) [57]. This requires 
that the weak scale values of bino and wino masses to be approximately equal. 
If instead these are roughly equal in magnitude but differ in sign, bino-wino 
mixing is suppressed, but agreement with the observed relic density is possible 
via bino-wino co-annihilation (BWCA) |58] . 

Of interest to us here is the potential for an enhanced rate for bottom quark 
production in SUSY events that occurs for MHDM, as exemplified by (but not lim- 
ited to) the HB /FP region of the mSUGRA model \59\ [5H] , or models where third 
generation squarks are significantly lighter than other squarks as, for instance, in 
the stop co-annihilation region of mSUGRA, in inverted hierarchy models where 
third generation sfermions are much lighter than those of the first two generations 
|60l [6T] , or in the framework suggested in Ref . [53j . 

It has been shown previously [38J that using 6-jet tagging techniques that are avail- 
able at the LHC, the SUSY reach may be enhanced by as much as 20% for parameters 
in the HB/FP of the mSUGRA model. Toward this end, we examine the reach of 
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the LHC with and without 6-jet tagging, in several models motivated by the rehc 
density measurement just discussed as well as by other considerations, to precisely 
delineate the circumstances under which 5-jet tagging will significantly enhance the 
LHC reach. Since SUSY events may also be enriched in t-jets, we also examine 
prospects for top jet tagging in SUSY events at the LHC. 

Having completed this introduction, for the benefit of the reader we provide a sum- 
marized version of the chapter's structure, which we have divided as follows: 



Sec 3.2 we introduce the various models chosen for the study. 



Sec 13.31 event simulation and calculational details are discussed. 



• Sec 3.4 we discuss the simulation of the signal and the analysis cuts. 



Sec 3.5 the results for the different models are presented. 



Sec 3.6 results specific to isolating third generation squarks are discussed. 



Sec 3.7 top-tagging is introduced, and its relative efficiency compared to b- 



tagging is discussed. 



Sec 3.8 charm tagging is discussed, including its relative merits in previous 



studies in comparison to the LHC. 



Sec 3.9 we summarize our work presented in this chapter. 



3.2 Models 

Here we discuss several models in which we may expect third generation fermions to 
be preferentially produced in SUSY models. We begin with the familiar mSUGRA 
model, and work our way through various other models motivated either by the relic 



density observation discussed in Sec. |3.1[ or by other considerations. 
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3.2.1 The mSUGRA model 

The mSUGRA model [62J was already introduced in Chap. ??. We remind the 
reader that this model is completely specified by the parameter set, 

tuq, mi/2, tan f3, Aq and sign(//) . (3-3) 



Typically, the weak scale value of |^| is similar in magnitude to rrig, and the bino is 
the LSP. However, for any chosen value of the requirement that electroweak 

symmetry be correctly broken imposes an upper bound on mo, since the value of 
/^■^ becomes negative for yet larger values of mo- There is thus a contour in the 
mo — w-i/2 plane where jj? = 0. For values of mo just below this bound, fj? <C m| 
and can be comparable to the SSB bino mass parameter. Mi so that the lightest 
neutralino is a mixed bino-higgsino state that can annihilate rapidly in the early 
universe, mainly via its higgsino content. This is the celebrated HB/FP region of 
the mSUGRA model [5lj, one of the regions of mSUGRA parameter space where 



the expected neutralino relic density is consistent with (3.2) |64] . For parameters in 
this region, squark masses are in the multi-TeV range, and the reach of the LHC is 
determined by final states from gluino pair production: although the higgsino-like 
chargino may be light, the mass difference m^ — m~^ is small so that leptons from 
its decays are too soft to increase the reach beyond that obtained via the E^^^^ 
signal from gluino pair production [65j. Since the LSP couples preferentially to the 
third family via its higgsino component, cascade decays of the gluino to third gen- 
eration fermions tend to be enhanced. As a result, the requirement of a 6-tagged 
jet in SUSY events reduces SM backgrounds and enhances the LHC reach by 15- 
20% beyond the reach via the inclusive Ej^^^^ channel in the HB/FP region of the 
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mSUGRA model [38]. We should also mention that the 6-jet multiplicity may also 
be enhanced in the mSUGRA model if third generation squarks happen to be light, 
either because of large bottom quark Yukawa couplings when tan/3 is large, or be- 
cause the At parameter happens to be just right so that m^^ <^ niq, and ti mainly 
decays via ti — )■ bWi and tZi, or ti — bWZi. 

3.2.2 Inverted mass hierarchy models 

The evidence for neutrino oscillations [66j and its interpretation in terms of neutrino 
masses provides strong motivation for considering SO (10) SUSY grand unified the- 
ories (GUTS) [67J. Each generation of matter (including the sterile neutrino) can 
be unified into a single 16 dimensional representation of SO{10) while the Higgs 
superfields Hu and are both contained in a single 10 dimensional representation, 
allowing for the unification of both gauge (and separately) Yukawa couplings. 
50(10) may either be directly broken to the SM gauge group, or by a two step pro- 
cess via an intermediate stage of SU (5) unification. The spontaneous breakdown of 
5*0(10) with the concomitant reduction of rank leaves an imprint on the SSB masses 
which is captured by one additional parameter with a weak scale magnitude 
but which can take either sign [68]. The model is then completely specified by the 
parameter set, 

mi6,mio,mi/2,M|,,tan/3,ylo and sign(/i) . (3.4) 

where we have assumed a common SSB mass parameter mig and a different param- 
eter niiQ for matter and Higgs fields in the 16 and 10 dimensional representations, 
respectively. The GUT scale SSB masses for MSSM fields then take the form [68], 
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2 2 2 




(3.5) 




u,d 



Unification of Yukawa couplings is possible for very large values of tan /3 j 69^ I70| . 
The SO(IO) framework that we have just introduced naturally allows a phenomeno- 
logically interesting class of models in which the matter sfermion mass order is 
inverted with the order for the corresponding fermions [60j. Specifically, in models 
with Yukawa coupling unification, the choice 



for the SSB parameters serves to drive third generation sfermion mass parameters 
to sub-TeV values, leaving first and second generation scalars as heavy as 2-3 TeV. 
A positive value of M|, < (mi6/3)^ is necessary to obtain radiative electroweak 
symmetry breaking [61] . The multi-TeV values of first and second generation scalar 
masses ameliorate the SUSY CP and flavour problems without destroying the SUSY 
resolution of the gauge hierarchy problem, since the fields with substantial direct 
couplings to the Higgs sector (gauginos and third generation scalars) have masses 
below the TeV scale. Because third generation sfermions are significantly lighter 
than their first/second generation cousins, we may expect that SUSY events are 
enriched in b- (and possibly t-) quark jets in this scenario. 

3.2.3 Non-Universal Higgs Mass Models 

Within the mSUGRA model, if niQ = m|^^(GUT) is smaller than or comparable 
to my2J ''^Hu ^^^s to a large negative value at the weak scale. The minimization 



Al = 2m?n = 4m?, 



16 



(3.6) 
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condition for the (tree level) Higgs scalar potential which reads 

mjr - mjr tan^ P Ml ..Ml , , 



(where the last approximation is valid for modest to large values of tan /3) , then im- 
plies that 1^1 > |-/Vfi^2| so that the LSP is essentially a bino, while the heavier -inos 
are mainly higgsino-like. A way of avoiding this conclusion is to choose ttu^^^ (GUT) 
to be so large that runs to small negative values at the weak scale. Within the 
mSUGRA model, this can only be realized by choosing mo ^ '^\i2 which gives us 
the well studied HB/FP region with MHDM discussed above. 

A different way would be to relax the assumed universality |50j between the matter 
scalar and Higgs boson SSB mass parameters in what has been dubbed as non- 
universal Higgs mass (NUHM) models, and adopt a large value for m|^^ (GUT). In 
order to avoid unwanted flavour changing neutral currents, we maintain a universal 
value mo for matter scalars. The GUT scale value of the SSB down Higgs mass pa- 
rameter may (may not) be equal to vn^H^ leading to a one (two) parameter extension 
of the mSUGRA framework that we wiU refer to as the NUHMl (NUHM2) model 
|51j . The NUHMl model is thus completely specified by the mSUGRA parameter 



set together with = sign{m jj^ J ^\mjj J, i.e. by, 

mo, m^, mi/2) ^0) tan/3 and sign(//) (NUHMl). (3-8) 

If is chosen to be sufficiently larger than mo, the parameter m^^ runs down 
to negative values but remains small in magnitude so that we obtain MHDM for 
any value of niQ and Curiously, the NUHMl model accommodates another 



possibility of getting agreement with (3.2). If m,^ < 0, fnjj^ and itt-h^ both run to 



^Of course, if m^ is chosen to be too large then rnj^^ does not run to negative 
values and electroweak symmetry breaking is no longer obtained. 
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large, negative values at the weak scale so that 

rnl = ^H^ + mjj^ + - mjj^ - mjj^ - M| (3.9) 

may be small enough for neutralinos to annihilate via the A and H resonances. 
Within the NUHMl framework, the Higgs funnel thus occurs for all values o/tan f3. 
Since the Higgs bosons A and H with relatively small masses are expected to be 
produced via cascade decays of gluinos and squarks, and since these decay prefer- 
entially to third generation fermions, we may once again expect an enhancement of 
the b- and, perhaps also, t-jet multiplicity. 

The NUHM2 model requires two more parameters than the mSUGRA framework 
for its complete specification. While these may be taken to be the GUT scale values 
of ni'j^^ and it is customary and more convenient to eliminate these in favour 

of niA and /x, and work with the hybrid parameter set, 

mo,mi/2,mA,/U, Ao,tan/3 (NUHM2) . (3.10) 

This then allows us to adjust the higgsino content of charginos and neutralinos at 
will, and furthermore allows as much freedom in the (tree-level) Higgs sector as in 
the unconstrained MSSM. 

3.2.4 Low iMal Dark Matter Model 

Instead of relaxing the universality between scalar masses as in the NUHM model, 
we can also relax the universality between the gaugino mass parameters. If we 
adjust the GUT scale value of M1/M2 so that Afi ~ M2 at the weak scale, we 
obtain mixed wino DM [57]. Since there is no principle that forces M1/M2 to be 
positive, we can instead adjust this ratio so that Mi ~ —M2 at the weak scale. In 
this case the LSP remains a bino with charged and neutral winos close in mass to it 
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and agreement with (3.2) is obtained via bino-wino co- annihilation |58| . Although 
collider signatures are indeed altered from mSUGRA expectations, we do not expect 
any enrichment of 6-jet multiplicity in this case. 



Although not obvious, agreement with (3.2 ) is also obtained if we maintain Mi = M2 



at Q = Mqut ) but instead reduce the value of | M3 1 . Specifically, for smaller values 
of IM3I, the (top)-squark mass parameters and also are driven to smaller values 
at the weak scale. These smaller values of top-squark masses and of A^, in turn, 
slow down the evolution of rrijj^ so that it runs to negative values more slowly than 
in the mSUGRA model. As a result, the weak scale value of rnjj^ though negative, 
has a smaller magnitude than in the mSUGRA case, so that the value of /x^ is 
correspondingly reduced [see Eq. ( |3.7| )] and the LSP becomes MHDM [55J. This is 
referred to as the low IM3I DM (LM3DM) model, and the corresponding parameter 
space is given by, 

mo,mi/2,M3, Ao,tan/3,sign(^) (LM3DM) . (3.11) 

Here mi/2 > denotes the GUT scale value of Mi = M2, while M3 (which is 
either positive or negative) denotes the corresponding value of M3 at the GUT 
scale. For niQ ~ rni/2 ~ 1 TeV, the GUT scale value of IM3I must be reduced 
from its mSUGRA value in order to obtain MHDM as discussed above. In contrast, 
if we fix mi/2 — 1 TeV, and take niQ to be multi-TeV, MHDM is obtained for 
values |M3|/mi/2 > 1- To simplify fine tuning issues, we will confine ourselves to 



mo ^ 1 TeV where we can obtain agreement with (3.2) by reducing the value of 
IM3I. We may expect an enhancement in the 6- multiplicity from SUSY events at 
the LHC because of the enhanced higgsino content of the LSP. 
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3.2.5 High M2 Dark Matter Model 

Very recently, it has been pointed out [56] that raising the GUT scale value of M2 
from its unified value of 772-1/2 to about (2.5-3)mi/2 for M2 > 0, or to between —2 
and —2.5 times mi/2 foi' -^2 < 0, also leads to a small value of |/i|, giving rise to a 
relic density in agreement with ( |3.2[ ). 

The parameter space of this high IM2I dark matter (HM2DM) model is given by, 

mo,mi/2,M2, Ao,tan/3,sign(/i) (HM2DM) . (3.12) 

where IM2I, the GUT scale value of the wino mass parameter, is dialled to large 
magnitudes to obtain MHDM. The large value of IM2I causes the Higgs SSB rnj^^ 
to initially increase from its GUT scale value of uiq as Q is reduced from Mgut- 
Ultimately, however, the usual top quark Yukawa coupling effects take over, causing 
m|^^ to evolve to negative values resulting in the well-known radiative breaking of 
electroweak symmetry. However, because of its initial upward evolution, the weak 
scale value of rn'j^^ is not as negative as in models with unified gaugino masses, and 
the value of fj? is correspondingly smaller. The neutralino LSP then has a significant 
higgsino component, and we may expect an enhancement of 6-jets in SUSY events 
at the LHC. 

3.3 Event simulation and calculational details 

We use IS A JET 7.74 [32] with the toy calorimeter described in Chap. [?] for the 
calculation of the SUSY signal as well as of SM backgrounds in the experimental 



environment of the LHC. We use parameters as described in Sec. 2.2.2 We con- 
servatively take the tagging efficiency ei, = 0.5 at the LHC design luminosity of 
100 fb^^/y, and assume that gluon and light quark jets can be rejected as b jets by 
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a factor R}, = 150 (50) if Et < 100 GeV {Et > 250 GeV) and a linear interpolation 
in between [SYj. For jets not tagged as a 6-jet, we require Exij) > 50 GeV. 

Gluino and squark production is the dominant sparticle production mecha- 
nism at the LHC for gluino and squark masses up to about 1.8 TeV, if rriq ~ rUg. If 
instead squarks are very heavy, gluino pair production will dominate the sparticle 
production rate up to about rrig ~ 0.8 TeV. Cascade decays of the parent gluinos 
and squarks then lead to signals in various multi-jet plus multi-lepton plus Ej^^^^ 
topologies [72] . 

Since SUSY particles are expected to be heavy (relative to SM particles) 
sparticle production is expected to be signalled by events with hard jets, possibly 
with hard, isolated leptons and large S™'^*^. The dominant physics backgrounds to 
these events with hard jets come from tt production, V + j production {V = W, Z), 
VV production and QCD production of light jets, where the Ej^^'^'^ comes from 
neutrinos produced by the decays of W or Z bosons or of heavy flavours. Missing 
Et may also arise from mismeasurement of jet or lepton transverse momenta and 
from uninstrumented regions of the detector. These non-physics sources of E^^^^^ are 
detector-dependent, and only qualitatively accounted for in our simulation with the 
toy calorimeter. With the hard cuts that we use to obtain the reach, we expect that 
the physics backgrounds will dominate the difiicult-to-simulate detector-dependent 
backgrounds, and the results of our analyses of the SUSY reach will be reliable. 
This expectation is indeed borne out since results of previous theoretical analyses 
of the SUSY reach [401 H2] compare well with the projected reaches obtained by the 
CMS [TU] and ATLAS [20] collaborations. The 51am in reach, if any, that we obtain 
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from 6-jet tagging, should if anything be more rehable than the absolute value of 
the reachJ3 

In the analysis detailed in the next chapter, we have examined the reach of the 
LHC for a wide range of sparticle masses, for the different models introduced in 



Sec. 3.2 To facilitate this, we generate signals and backgrounds (calculational details 
are described below) and only write out events that include at least two jets with 
EtU) > 100 GeV and E^^^^ > 100 GeV, which we refer to as our basic cuts. 
The corresponding cross sections for SM events are shown in the second column of 
Table |3.1[ For low to medium values of sparticle masses, the sparticle production 
cross sections are large enough for us to extract the signal above SM backgrounds 
with relatively soft analysis cuts. For very heavy sparticles, however, the production 
rate is small, but essentially all events contain very energetic jets and large i?™^®. 
The detection of the signal is then optimized by using very hard cuts that strongly 
suppress SM backgrounds while retaining bulk of the SUSY signal. Since our aim is 
to develop a strategy that can be applied to essentially the entire interesting mass 
range of a wide variety of models, we are led to evaluate the signal together with 
the SM background for a wide range of cuts, detailed in the next section. 
To understand the relative importance of the different background sources, in the 
last three columns of Table [ST] we list the corresponding cross sections for the softest 



set of cuts that we use in our analysis detailed in Sec. |3.4 

In the last two rows we also list the corresponding signal cross sections for two 
WMAP-consistent cases in the HB/FP region of the mSUGRA model. Several 



^The absolute reach may also suffer from the fact that SM backgrounds may 
be somewhat larger than those obtained using shower Monte-Carlo programs when 
proper matrix elements for jet production are included. We expect though that the 
gain in the reach from 6-tagging may again be less sensitive to the inclusion of the 
proper matrix elements. 
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iSonrrp 


'^basic 


^cut y^^ ) 


^ cnx,\^^ } 


^cut V y 


ff 


IQQDD 

X U \J\J\J 


2.16 


1.41 


365 


W + j 


21400 


12.0 


1.36 


0.133 


z + j 


8850 


5.11 


0.059 


0.0052 


VV 


89.8 


0.0248 


0.0020 


0.0001 


QCD 


93700 


11.6 


3.11 


0.467 


Total 


1.44 X 10^ 


30.9 


5.94 


0.97 


mSUGRAl 


261 


12.0 


9.26 


3.86 


mSUGRA2 


48.4 


2.44 


1.95 


0.87 



Table 3.1: Cross sections in fb for the SM production of tt, W + j, Z + j, VV, and 
QCD jet events that form the dominant backgrounds. The second column gives the 
cross section for events with the basic requirements of two jets with ExiJ) > 100 GeV 
and Ej^^^^ > 100 GeV. The last three columns give the corresponding cross sections 



for the softest set of cuts listed in Table 3.2 and Table [373] that we actually use in 
our analysis, with no requirement of b-jet tagging (column 3), requiring at least one 
tagged 6-jet (column 4) and at least two tagged 6-jets (column 5). For illustration, we 
also list the corresponding signal cross sections for two points in the HB/FP region 
of the mSUGRA model, with Aq = 0, tan /3 = 10 and ~ 1 TeV, rUg ~ 3 TeV 
(mSUGRAl) and rtig ~ 1.5 TeV and rrig ~ 3.9 TeV (mSUGRA2). 



comments are worth noting. 



• We see that with the basic requirements of two jets with Et > 100 GeV and 
j^miss > ]^QQ QeV, the background is two (three) orders of magnitude larger 
than the signal for nig ~ 1 (1.5) TeV; however, the analysis cuts very efficiently 
reduce the background, while reducing the signal by a much smaller factor. 

• After these analysis cuts we see that QCD, followed hy V + j production, are 
the leading backgrounds to the inclusive E^^^^ signal. Top pair production, 
while significant, is considerably smaller. Since we do not require the presence 
of leptons, the background from VV production is negligible. 

• The backgrounds from QCD and V + j production can be sharply reduced by 
the use of 6-jet tagging with relatively small loss of the signal. In contrast. 
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since top events necessarily contain 6-jets, 6-tagging reduces the tt background 
only by a modest amount. 



Table 3.1 highlights the importance of a careful evaluation of the QCD and 
the V + j backgrounds. This is technically complicated because the large size of 
the cross sections necessitates simulations of very large number of events to obtain 
a reliable estimate for the backgrounds after the very hard cuts that are needed 
for optimizing the reach of the LHCj^ Moreover, since the cross section is a rapidly 
falling function of the centre of mass energy, or equivalently, the hard scattering pT 
of the initial partons, we must ensure that our procedure generates events even for 
very large values of where the matrix element is very small, resulting in a much 
smaller weight. To facilitate this, we have generated the various backgrounds using 
different numbers, Nl^^ , of hard scattering bins: the bin intervals are finely spaced 
for low values of where event weights are very large. We choose 7V/^^ = 53, 13, 8 
and 7 for i = QCD, V + j, tt and VV, respectively, where the choice N^^^ = 53 
reflects the largeness of the QCD cross section. We have generated a total of about 
lOM QCD events, about IM W + j events and about 500K-700K events for each of 
the other backgrounds. If, for any set of cuts, we find zero events in our simulation 
of a particular background, we set this background cross section to a value corre- 
sponding to the one event level in the bin with the smallest weight in our simulation. 



■^Of course, the fact that we are far into the tails of these backgrounds where 
the simulations (which will be tuned to data when these become available) require 
possibly unjustified extrapolations is a different matter. 
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3.4 Bottom jet tagging and the reach of the LHC 
3.4.1 Simulation of the Signal and the LHC reach 

Simulation of the signal events is technically much easier than that of the back- 
ground. This is largely because the signal typically originates in heavy sparticles, 
and so passes the hard analysis cuts with relative ease compared to the background. 
To assess how much 5-jet tagging extends the SUSY reach of any particular model, 
rather than perform extensive and time-consuming scans of the parameter space, we 
have defined "model lines" along which the sparticle mass scale increases. We then 
choose parameters along these lines, and for every such parameter set use IS A JET 
7.74 to generate a SUSY event sample. Next, we pass this event sample through 
the set of analysis cuts defined below, and define the signal to be observable at the 
LHC if for any choice of cuts 

• the signal exceeds 10 events, assuming an integrated luminosity of 100 fb^^, 

• the statistical significance of the signal A^signai/\/-^back > 5, and 

• the signal to background ratio, A'signaiZ-^back ^ 0.25. 

We also require a minimum of 15 events after cuts in our simulation of the signal. 
We obtain the reach for each model line by comparing the corresponding signal with 
the background, and ascertaining where the signal just fails our observability criteria 



for the entire set of cuts in Table and Table 3.3 An important part of our work 
involves tagging of b-jets at the LHC, using both the presence or complete absence 
of these tagged jets to achieve our results. 
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3.4.2 Analysis cuts 

The inverted mass hierarchy model based on 5*0(10) SUSY GUTs, whose hahmark is 
the hght third generation, serves as the prototypical case where we expect enhanced 
6-jet multiplicity in SUSY events. We have used this framework to guide us to the 
set of analysis cuts that can be used for the optimization of the SUSY signal for 
a wide range of sparticle masses in a wide class of models. Toward this end, we 
fix n < 0, Aq < 0, and tan/3 = 47 (a large value is needed for the unification of 
Yukawa couplings) and choose mio = V^mie, Aq = —2mie to obtain the hierarchy 
between the first/second and third generation scalars as discussed above. The choice 
Md = 0.25mi6 facilitates electroweak symmetry breaking. We vary the gluino mass 
along the "model line" with mi/2 = 0.36mo + 48 GeV which maintains a hierarchy 
between the generations. 
The value of 




is typically around 3.5-4.1 along this model line. 

The optimal choice of cuts depends on the (a priori unknown) sparticle spectrum, 
and to a smaller extent on their decay patterns. While hard cuts optimize the signal 
if sparticles are heavy, these would drastically reduce (or even eliminate) the signal 
if sparticles happen to be light. In order to obtain a general strategy that can be 
used for a wide variety of models, we have used the SO{10) model with < to 
define a universal set of cuts that can be used for SUSY discovery in any of the 
various models that we have introduced, and likely, also for a wider class of models. 
Toward this end, we generate a sample of signal events for this "test model line" 
and run this, as well as the SM backgrounds that we discussed above, through each 
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one of the large set {i.e. the complete set that includes the numbers listed in the 



parentheses) of analysis cuts detailed in the nine rows in both Table 3.2 and Ta- 



ble 3.3 Here, rn-gg is the scalar sum of the transverse energies of the four hardest 
jets in the event combined with the missing transverse energy, A(j) is the transverse 
plane opening angle between the two hardest jets, and A(j)b the corresponding an- 
gle between the two tagged 6-jets in events with rih > 2. To clarify, the softest 
set of cuts that we use for the 06 signal has [i?™^^, EtUi), £'t(J2); -E't(6i)i ^es] 
> (300,300,200,40, 1500) GeV, nj > 4 and transverse sphericity St > 0.1, with no 
restriction on jet opening angles. Note that because Et^Ji) > ET{j2), there are 21 
combinations for the minimum values of [-©^'(ji), -Et(J2)] that we have used. Next, 
we harden the cut on one of these observables to the next level, keeping the others 
at the same value, etc. until the complete set of 6x21x5x6x3x3x5x2 com- 
binations has been examined for > 2. Since there are (is) no (just one) tagged b 
jets in the nj, = (ub = 1) case, there are correspondingly fewer combinations for 
these analyses. 

For each of these cut choices, we analysed the observability and statistical 
significance of the LHC signal for our test 5*0(10) model line for an integrated lu- 
minosity of 100 fb~^. We found that a subset of cuts was sufficient to ensure the 
observability of the SUSY signal over the entire mass range. Specifically, restricting 
the minimum values of the transverse energies of the two hardest jets to the eleven 



combinations shown in Table 3.4, and dropping the cut values shown in parenthesis 



for the other variables in both Table |3.2| and Table 3.3 had no impact upon the 
observability (and the statistical significance) of the signal over the entire sparticle 
mass range. 
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Variable 06, 16 

Ef'^ (GeV) > (300), 450, ...,900, (1050) 

-Er(ji) (GeV) > 300, 400, 800 

-Et(J2) (GeV) > 200, 300, Ef''{ji) - 100 GeV 

ETih) (GeV) > 40, 100, 200, (300, 400) 

meff (GeV) > 1500, 2000, 2500, 4000 

Acj) < 180°, 160°, 140° 

A^b < N/A 

Hj > 4, 5, 8 

St > 0.1, (0.2) 



Tabic 3.2: The complete set of cuts examined for extraction of the SUSY signal over 
the SM backgrounds. The 06 and 16 entries denote requirements for events without 
any restriction 6-jet tagging, or with at least one tagged 6-jet. 



Variable 26 

Ef''' (GeV) > (300), 450, 600, 750, (900, 1050) 

EtUi) (GeV) > 300, 400, 800 

^t(J2) (GeV) > 200, 300, Ef'^'iji) - 100 GeV 

ETih) (GeV) > 40,100,200,300,(400) 

meff (GeV) > 1500, 1750, 2000, 2250, (2500, 2750) 

Ac/) < 180°, (160°, 140°) 

A(j)b < 180°, 150°, 120° 

nj> 4, 5,. ..,7, (8) 

St > 0.1, (0.2) 



Table 3.3: The complete set of cuts examined for extraction of the SUSY signal 

over the SM backgrounds. The 26 entries denote requirements for events with at 
least two tagged 6-jets. For the final analysis of the reach in the various models, we 
dropped the cut values within the parenthesis, and replaced the 21 combinations 
for the minimum values of [ET{ji), ET{j2)] with the 11 combinations in the last two 
rows (below the horizontal lines) in the table. 
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[ET{ji),ET{j2)] (GeV) > (300, 200), (400, 200), (500, 200), (500, 300), 

(500, 400), (600, 200), (600, 500), (700, 300), 
(700, 600), (800, 300), (800, 600) 



Table 3.4: For the final analysis of the reach in the various models, we dropped the 
cut values within the parenthesis in Table [3^2] and Table 3.3 and replaced the 21 
combinations for the minimum values of [ETiji), ET{j2)] with the 11 combinations 
in the table above. 



In the remainder of this project we, therefore, confine ourselves to this lim- 
ited subset of cuts, as this speeds up the analysis considerably. 



3.5 Results for LHC reach using 6-jet tagging 

In this section, we evaluate prospects for increasing the reach of the LHC by the 
use of 6-tagging to reduce SM backgrounds, thereby increasing the statistical signif- 



icance of the SUSY signal, for each of the models introduced in Sec. 3.2 We confine 
ourselves to various 1-parameter model lines (introduced below) along which sparti- 
cle masses increase and run the signal and backgrounds through each of the final set 



of cuts in Table 3.2 and Table 3.3 and optimize the signal by selecting the cut choice 
that yields an observable signal with the highest statistical significance. To assess 
the gain from 6-tagging, for each model line we first do so without any requirement 
on 6-tagging, and then repeat it requiring, in addition, at least one and at least two 
tagged 6-jets. 
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3.5.1 The HB/FP region of the mSUGRA model 

The possibility of increasing the LHC reach was first studied in the HB/FP region of 
the mSUGRA framework [38] , where it was found that the reach could be increased 
by up to 15-20%. We have repeated this study, albeit with a somewhat different 
model line with 

mi/2 = 0.295mo - 507.5 GeV,tan/3 = 30, = 0, 



in the HB/FP region that saturates the relic density in (3.2) and of course, with the 
different set of cuts that we use here. We find an increased reach from 6-tagging in 
qualitative agreement with Ref. [38]. 

3.5.2 Inverted mass hierarchy model 



As discussed in Sec. 3.4.2 we have already used the 50(10) model with /i < and 



parameters related by (3.6) where we obtain an inverted mass hierarchy to choose 
the final set of cuts for our analysis. Here, we show results for the reach of the LHC 
with and without requirements of 6-jet tagging for two model lines with a significant 
inversion of the sfermion mass hierarchy, one for each sign of /U. For both of these, 
we choose 

- Ao = 2mi6 = \/2mio, tan /3 = 47, (3.13) 



with 



Md = 0.25mi6 and mi/2 = 0.36mi6 + 48 GeV for n < 0, (3.14) 
Md = 0.20mi6 and mi/2 = 0.30mi6 + 39 GeV for ^ > 0. (3.15) 



Our results are shown in Fig. 3.1, where we plot the largest statistical sig- 



nificance of the signal, A'signai/ V ^back, versus the corresponding gluino mass for 
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(a) /i < 0, and (&) ^ > 0, assuming an integrated luminosity of 100 fb ^. The maxi 
mal A'signai/ V -^back was obtained running over all the cuts in Table 



3.2 



and Table 



3.3 



subjecting it to the requirement that the A'signaiZ-^back > 0.25 and A'^signai > 10 event 
criteria are satisfied. The solid (red) curves show this significance for the inclusive 
Ej^^^^ signal with no requirement of 6-jet tagging, while the dashed (black) curve 
and the dotted (blue) curves correspond to cases where we require at least one and 
two tagged 6-jets, respectively. The wiggles in these curves reflect the statistical 
errors in our simulation. 

We attribute the somewhat larger reach in the left frame to the fact that the mass 
hierarchy (as measured by the value of S) is somewhat smaller for ^ < 0, so that 
qg makes a larger contribution in this case. We also see that for fi < 0, 6-tagging 
leads to an increase of the LHC reach by ~ 200 GeV, or about 10%, while the 
corresponding increase is somewhat smaller for the model line with positive /i. 
This difference (which may well not be very significant in view of the wiggles) is 
evidently due to the increased reach in the 2b channel, and could arise from a com- 
plicated interplay between the effect of cuts and the sparticle spectrum: for instance, 
for rrig ~ 1960 GeV, nij^^ is significantly lighter in the ;U < case, while m^^ is con- 
siderably heavier. As a result, the branching fraction for the decays g ^ bbi, which 
likely leads to a harder spectrum for 6-jets (compared to 5 — )■ tti, which constitutes 
the bulk of the remaining decays of the gluino), falls from 38% for negative to 
28% for positive /i. 
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Figure 3.1: The statistical significance for the inverted hierarchy SO{W) model hnes 
introduced in the text, assuming an integrated luminosity of 100 fb^^ for (a) /i < 0, 
and (b) ^ > 0. The sohd (red) hne is for the signal with no requirement on 6-tagging, 
the dashed (black) line is with the requirement of at least one tagged 6-jet, and the 
dotted (blue) line is with at least two tagged 6-jets. The signal is observable if the 
statistical significance is above the horizontal line at A^'signai/V-^back = 5. 

Non-universal Higgs mass models 

Next, we turn to the impact of 6-tagging on the reach in NUHM models with 
just one additional parameter that is adjusted so that agreement with the ob- 
served relic density is obtained either by tempering the LSP content so that it is 
MHDM (m^ > mo), or by adjusting the masses so that the LSP annihilation rate 
is resonantly enhanced by the exchange of neutral A or H bosons in the s-channel 
{■171(1) < 0). We did not study the NUHM model where both Higgs SSB mass param- 
eters are arbitrary - the so-called NUHM2 models in the nomenclature of Ref. |51| 
- because this meant that both ttia and fi are arbitrary, resulting in too much free- 
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dom for definitive analysis. Beginning with the MHDM cases of the LSP where 
sparticle decays to third generation quarks are enhanced by the higgsino content 
of the LSP, we introduce two model lines with Aq = 0, tan/3 = 10 and /U > 0, 
with (1) niQ = mi/2, and (2) niQ = 3mi/2- In the former case, the squarks of the 
first two generations are roughly degenerate with gluinos, whereas in the latter case 
rriq ~ 1.6mg. 



Our results for the statistical significance of the LHC SUSY signal, with and 



without 6-jet tagging are shown in Fig. 3.2 for (a) rriQ = mi/2, and (b) mo = 3mi/2- 
We see that while 6-tagging clearly improves the reach by ~ 10% in the case shown 
in frame (b), it leads to a degradation of the reach in frame ( a) . 
We have traced this to the fact that for this case where squark and gluino masses 
are comparable, squark production (particularly first generation squark production) 
makes a significant contribution to the signal after the hard cuts. Then, since un- 
like gluinos which decay "democratically", these squarks decay to charginos and 
neutralinos (remember that because rriq ~ nig, the decay q ^ qg is suppressed by 
phase-space) plus quarks of their own generation, a sizeable fraction of the inclusive 
S™'^^ signal is actually cut out by any 6-tagging requirement. 

In frame ( b) , the squarks are much heavier than gluinos and so contribute a smaller 
fraction of the signal, but more relevantly, q ^ qg with a large branching fraction, so 
that 6-tagging helps in this case. These considerations also explain why the increase 
in reach from 6-tagging is not as large as in the case of the HB/FP region of the 
mSUGRA model where mg ^ rUg [38] . 



We now turn to the < model line shown in Fig. 3.2 ; for which we have chosen 
ruQ = 5mi/2 (to ensure squark contributions to the signal do not dilute the effect 
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of b tagging as in the case that we just discussed), Aq = 0, tan/3 = 20 and fJ, > 0, 



and < is adjusted to give agreement with (3.2) via resonant annihilation of 
LSPs through A/H exchanges in the s-channel. This means that A and H must be 
relatively light and accessible in cascade decays of gluinos and squarks. However, 
we see no enhancement of the LHC reach in this case. 

We understand this in hindsight. In this case is large so the lighter neutrali- 
nos produced in gluino cascade decays are gaugino-like, with rn^^ ~ ~ ^"^Zi" 
Then the very condition 2m~ ~ tjia that makes the LSP anihilation cross section 
resonant suppresses the phase space for the decays of Z2 ^ A or H + Zi, so that 
these are not significantly produced in cascade decays of gluinos. Since squarks are 
very heavy, they are essentially irrelevant to this discussion. 



Low M3 dark matter model 

As explained above, we can also obtain MHDM, and hence a potential increase in 
reach via 6-tagging, in models with non-universal gaugino mass parameters where 
|-/Vf3(GUT)| is taken to be reduced compared to its value in models with gaugino 
mass unification. To study the gain in the reach that we may obtain in this case, 
we have explored an LM3DM model line with 

mo = mi/2,Ao = 0, tan/3 = 10, ^ > 0, 

where the GUT scale value of M3 (which we take to be positive) is adjusted to 
saturate the measured CDM relic density]^ The corresponding dependence of the 



statistical significance of the SUSY signal on mg is shown in Fig. 3.3 We see that 



"^Roughly speaking, for mo = m^/a = 700 GeV, M3(GUT) = 277 GeV, and for 
an increase of 6mQ in mo = mi/2) the GUT scale value of M3 has to be raised by 
about 6Ms ~ 5mo/2.25. 
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in this case 6-tagging leads to an increase in reach close to 15%. This is because 
though gluinos and squarks are both reduced in mass relative to their uncoloured 
cousins, the reduced value of the gluino mass parameter leads to rUg ~ (1.4 — 1.5)mg 
even for mo = to be compared to niq ~ rrig that we obtained for models with 

unified gaugino masses as e.g. in the NUHM case just discussed. The large value of 
rriq relative to rrig then leads to an enhanced reach via 6-tagging just as before. 

High M2 dark matter model 

As a final example, we consider the LHC reach in the HM2DM model, where agree- 
ment with (3.2) is obtained by raising |M2(GUT)| from its canonical value of 



in models with gaugino mass unification, so that the lightest neutralino is MHDM. 
Since the LSP contains a substantial higgsino component, it is again reasonable to 
expect that 6-jet tagging may increase the SUSY reach of the LHC. 
As we have already seen in other examples, the increased reach from 6-jet tagging 
depends on the value of the squark mass relative to rUg. This led us to consider 
two model lines with, (a) uiq = mi/2, and (6) uiq = M2(GUT), for both of which 
we take tan/3 = 10, /i > and Ao = 0. Since the correct relic density is obtained 
by raising M2, model-line (b) which gives heavier squarks than model line (a) will 
give a smaller reach as measured in terms of rrig. The increase in the reach from 
6-jet tagging will, however, be larger for model line (b) since squark contributions 
to sparticle production are kinematically suppressed. 

The statistical significance of the SUSY signal in the HM2DM model is shown 



for the two model lines in the two frames of Fig. 3.4 Indeed we see that while the 
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reach in the left frame for uiq = mi/2 extends to rrig < 2.5 TeV (as compared to 
2.1 TeV in the right frame), there is very httle gain in the reach from 6-jet tagging 
in this case where squark and gluino masses are comparable. This is in contrast to 
the gain in reach of ~ 8% for the case of heavier squarks in the right hand frame. 

3.6 Is direct detection of third generation squarks pos- 
sible? 

Establishing that any new physics signals at the LHC arise from supersymmetry 
will require the identification of several superpartners. In models where the third 
generation is significantly lighter than the other generations, it is natural to ask 
whether it is possible to detect signals from the direct production of third generation 
squarks. As already mentioned, their detection as secondaries from production and 
subsequent decays of gluinos is possible if the gluino itself is not very heavy [73] . Our 
goal, therefore, is to examine whether the signal from the direct production of third 
generation squarks can be separated both from SM backgrounds, as well as from 
production of other SUSY particles. Clearly, this is a model-dependent question, 
since the SUSY "contamination" to the third generation signal will depend strongly 
on the masses of the other squarks and the gluino. In this section, we will study 
this issue within the context of the inverted mass hierarchy model with /i < 0, that 
we have used as our canonical test case. 

Since there are essentially no third generation quarks in the proton, the cross 
section for third generation squarks falls rapidly with the squark mass, and the 
signal becomes rapidly rate-limited. We show this with the solid (blue) curve in 
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Fig. 3.5 where we plot the cross section in fb vs. the squark mass in GeV units, and 
contrast the results for the third generation squarks with those from either all gen- 
eration squarks only, or from gluinos plus all generation squarks for niq = nig. Here 
we show the cross sections corresponding to production of all generation squarks as 
the dotted (purple) curve, and all generation squarks and gluinos as the dashed- 
dotted (black) curve. We see that the cross section for third generation squarks 
is a subdominant part of the SUSY cross section, and further that it is very low, 
approaching 1 fb. 

Therefore, we confine ourselves to the signal from third generation squarks with 
masses around 300-500 GeV, where the signal is likely to be the largest. To un- 
equivocally separate out the third generation signal, we must use cuts that are hard 
enough to reduce the SM backgrounds to acceptable levels, yet not so hard as to en- 
hance the "contamination" from heavier sparticles that, though they are produced 
with (much) smaller cross sections than third generation squarks, would pass these 
hard cuts with much larger efficiency. 

Since third generation sfermions decay preferentially to third generation fermions 
(we focus on the case where ti — )• bWi is accessible), we study the signal with at 
least one tagged 6-jet. We found, however, that even the softest set of cuts in both 



Table 3.2 and Table 3.3 that we actually use for our analysis of the SUSY 6-tagged 
signal, are too hard for the purpose of extracting the signal from third generation 
squarks. 

We, therefore, returned to our basic cuts, 

Ef''''' > 100 GeY,ET{ji,j2) > 100 GeV 
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and augmented these with the requirements, 

ET{h,k) > 100 GeV, St > 0.1, > 1, 



and ran the third generation signal through the analysis cuts in Table 3.5 to extract 
the optimal A'signaiZ-^back ratio (where the background includes the SM and the 
SUSY contamination as we discussed). These cuts, which are applied "from below", 
primarily serve to control the SM background which is very large after just the basic 



cuts (see Table 3.1), but reduced by the additional requirements of a tagged b-jet 



and two additional 100 GeV jets. 



Variable Values 
nb> 1 
Ef''''' (GeV) > 100, 150, 200, 250 

[ET{ji),ET{j2)] (GeV) > (100, 100), (200, 100), (200, 150) 

(300, 100), (300, 150), (300, 200), 
(400, 100), (400, 150), (400, 200) 
Exih) (GeV) > 40, 100, 200, 300, 400 

rries (GeV) > 500, 600, 1500 

Uj > 4, 5, 6, 7 



Table 3.5: The set of cuts examined for the extraction of the third generation squark 
signal at the LHC. See the text for the additional cuts we have imposed. 



We show the results of our analysis in Table 3.6 The parameters are shown 
in the first four rows of the Table, while the next few rows show representative 
sparticle masses. The first two cases are along the /i < model line that we had 
introduced previously. In the first two cases B(ti — ?■ bWi) = 1, while in Case 3, 
-6(^1 —7- bWi) = 0.74, with the remainder being made up by the decay ti — )■ tZi. 
The next several rows list the optimized choice of cuts from the 4x9x5x11x4 



possibilities in Table 3.5, along with the cross sections for (i) the third generation 
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signal, (ii) the SM background, and {in) the "SUSY contamination" defined as the 
SUSY signal from production of sparticles other than third generation squarks, after 
these cuts. 

We see from these cross sections that both the event rates and the statistical signif- 
icance of the third generation signal (even with the SUSY contamination included 
in the background) is very large. The problem, however, is that the signal to back- 
ground ratio is smaller than 0.1, if the SUSY contamination is included in the back- 
ground, and fails to satisfy our observability criterion]^ We can, however, reduce the 
SUSY contamination (primarily from heavier sparticles) relative to the third gener- 
ation signal by requiring that the signal is not too hard. Toward this end, we impose 
an upper limit, rUcs < 1000 GeV, which efficiently reduces the contamination from 
heavy sparticles with correspondingly modest reduction of the cross sections from 
the softer third generation and SM processes. The corresponding cross sections after 
this cut are shown on the next three rows of the Table, while the last row shows the 
final two signal to total background ratio that we are able to obtain, along with the 
statistical significance of the third generation signal with an integrated luminosity 
of 100 fb~^ 

^Many authors do not impose such a requirement on the observability of the 
signal. We believe that some requirement on the A'signaiZ-^back ratio is necessary 
since otherwise a signal with 5K events, above a background of IM events would 
be considered significant. This would be indeed be the case if the background were 
known to a very high precision; however a systematic uncertainty of 0.5% on the 
background could clearly wipe out the signal, at least if the signal is extracted by 
subtracting the theoretically calculated background! In the case at hand, where the 
SUSY model is not a priori known, and has to be arrived at using the same data, 
it is clear that subtraction of the SUSY contamination will suffer from consider- 
able uncertainty until the data and theory both become mature enough for such a 
subtraction to be carried out. While our criterion requiring A'signaiZ-^back > 0.25 is 
admittedly arbitrary, we believe that it is necessary to impose some lower limit on 
the signal to background ratio for a semi-realistic assessment. 
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Several comments about the Table are worth noting. 

• We see from the Table that before the cut restricting the value of niee from 
above, the background was dominated by SUSY contamination. In contrast, 
after this cut, the dominant source to the background comes from SM pro- 
cesses. 

• With the cuts that we have devised, the event rates for the third genera- 
tion signal as well as its statistical significance are large. For reasons already 
discussed, we do not, however, believe that it will be easy to unequivocally 
ascertain the direct production of third generation squarks in the signal. For 
this to be unambiguously possible, it will be necessary to have an understand- 
ing of the contributions from other SUSY sources to the event rate after our 
cuts. This may well be possible because with hard cuts it should be possi- 
ble to isolate the signal from heavy squarks and gluinos where contamination 
from both SM and the lighter third generation squarks is small. Just how well 
it will be possible to extrapolate this measured signal into "softer kinematic 
regions" will determine the precision with which the SUSY contamination can 
be subtracted. This issue is beyond the scope of the present analysis 

• We examined additional cuts on ET{ji,j2) and rij to see if we could raise 
the signal to background ratio. We found that a small increase 10%) may 
indeed be possible by restricting nj from above to be smaller than 8 or 9. 
Since our calculation of the background with high jet multiplicity is carried 
out only in the shower approximation, we did not feel that our estimate of this 
improvement is reliable, and choose not to include it in the Table. 
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• We stress again that the SUSY contamination is model-dependent. We can 
see from the Table that if gluinos and other squarks are indeed decoupled at 
the LHC, and only third generation squarks are light, their signal should be 
readily observable in all three cases. 

3.7 Top tagging and the reach of the LHC 

We have seen that requiring a 6-tagged jet reduces the SM background relative to 
the SUSY signal in a wide variety of models, and so increases the SUSY reach of 
the LHC. This then raises the question whether it is possible to further increase 
this reach by requiring a top-tagged jet, since the mechanisms that serve to enhance 
the decays of SUSY particles to 6-quarks frequently tend to enhance decays to the 
entire third generation. SM backgrounds to E^^^^ events with t-quarks should, of 
course, be smaller than those for events with 6-quarks. In this section, we study the 
prospects for top tagging, once again using the inverted mass hierarchy model line 



(3.14) to guide our thinking. 



Top tagging in SUSY events has been suggested previously for the reconstruc- 
tion of SUSY events, assuming that ti or bi are light enough so that g — )• tti — )• tbWi 
and/or bbi — ?• btWi occur with large branching fractions [73]. It was shown that for 
rrig ~ 700 GeV, for which the SUSY event rate is very large, partial reconstruction 
of SUSY events with gluinos decaying to third generation squarks was possible at 
the LHC. 

We follow the approach developed in this study to reconstruct the top quark 
via its hadronic decay mode. In a sample of multi-jet + events with at 
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least one tagged 6-jet, we identified a hadronically decaying top by first identifying 



all pairs of jets (constructed from those jets that are not tagged as a 5-jets) as a 
hadronically decaying W if \mjj — Mw\ < 15 GeV. We then pair each such W with 
the tagged 6-jet(s) and identify any combination as a top if \mb\y — mt\ < 30 GeV. 
If we can reconstruct such a "top", we defined the event to be a top-tagged event. 
The efficiency for tagging tops in this way turns out to be smallj^ 
For our examination of the impact of top tagging on the SUSY reach of the LHC, 



we have chosen the SO{10) model line (3.14) with /i < as a test case. In this 
case, since other squarks are heavy, the gluino mainly decays with roughly equal 
likelihood via g — )• tit and g — )■ bib, where subsequent decays of the third generation 
squarks can lead to yet more top quarks in SUSY events. As for the case of b-jet 



tagging, we have run the SUSY sample through a set of cuts shown in Table 3.7 



to optimize our top-tagged signal relative to SM background. Because of the small 
efficiency for top-tagging we cannot, however, afford a large reduction of the signal 
from multiple cuts. We have, therefore, restricted our optimization to cuts on just 
the three variables Ej^^^^ , and nj , imposing the basic requirements on the signal 



^In a simulated sample of about 90K tt pairs with a hard scattering Et between 
50-400 GeV, we found only 6,255 top tags even with ei, = 1. To understand this large 
loss of efficiency we note that first, leptonically decaying tops (branching fraction 
of ~ 1/3) are clearly not identified. Second, 6-jets are within their fiducial region 
{Etj > 40 GeV, \r]j\ < 1.5, with a 5-hadron with pt{B) > 15 GeV within a cone of 
Ai? = 0.5 of the jet axis) only about 5/8 of the time. Third, it is necessary for the 
top with the 6-jet inside the fiducial region to decay hadronically in order to make 
the top mass window, since the wrong combination mostly falls outside. Finally, if 
the jets from the W from the top with the tagged b merge or radiate a separate jet 
at a large angle, this W is lost, and hence the top, is not tagged. We have checked 
with our synthetic top sample that the choice of mass bins of ±15 GeV about M\y 
and ±30 GeV about mt suggested in Ref. |73j does not lead to loss of signal from 
events where the top decays hadronically into well separated jets: most of the loss 
in efficiency comes from the other factors detailed above. 
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as discussed in Sec. 3^ The results of our SUSY reach analysis with top-tagging 
are summarized in Table 13.81 

Here, we show the optimized statistical significance of the SUSY signal for three 
cases in the vicinity of the ultimate reach using this technique. In this table, we show 
representative sparticle masses along with branching fractions for sparticle decays 
that lead to top quark production in SUSY cascades. We then detail the final choice 
of cuts that optimizes the top-tagged SUSY signal. We also show the top-tagged 
signal cross section after these cuts along with the corresponding SM background, 
and the statistical significance of the top-tagged signal achieved in cases 1 and 2; 
for case 3, the signal is not observable by our criteria. Finally, in the last two rows 
we show the corresponding statistical significance using 6-jet tagging discussed in 



Sec. 3.4, We see from the Table that while top tagging allows an LHC reach for 
nig just above 1600 GeV, the top-tagged rate becomes too low for heavier gluinos. In 
contrast, 6-jet tagging yields a statistical significance in excess of 50 close to the top- 
tagged reach. We thus conclude that while top-tagging can be used as a diagnostic 
tool, or even for reconstruction of SUSY events |73] in favourable cases, it will not 
extend the SUSY reach of the LHC. 



3.8 Charm-jet tagging 

Charm jet tagging offers a different possibility for enhancing the SUSY signal, espe- 
cially in the case where a light top squark dominantly decays via ti — t- cZi. Charm 
jets may be tagged via the detection of a soft muon within the jet. Muons inside 
jets also arise from semi-leptonic decays of 6-quarks and from accidental overlaps of 
unrelated muons with jets. Since is significantly larger than rric, the variables 
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Ip-p^l = Ipt(^) X Pj\ and Ai?(/u,j) = A(j){fi, j)"^ + A?7(/z, j)^ can serve to distin- 
guish muon-tagged c-jets from correspondingly tagged 6-jets or accidental overlap 
of an unrelated muon with jets. 

Charm jet tagging with soft muons was first examined in Ref. [74j as a way of en- 
hancing the i-squark signal from pp — )• ttX cc + Ej^^^^ + X production at Run 
I of the Fermilab Tevatron, but was found to have a reach smaller than the reach 
obtained via the conventional E^^^^ analysis because the muon-tagged signal was 
severely rate-limited. It was, however, subsequently shown that using soft muons to 
tag the c-jet indeed enhances the top squark reach [75] but only for an integrated lu- 
minosity larger than ~ 1 fb^^, available today after the upgrade of the Main Injector. 

These considerations led us to examine whether charm tagging may be sim- 
ilarly used at the LHC, at least for the case where t — )• cZi. Since the goal is to 
separate the charm jets from the decay of ti from other SUSY sources (which are fre- 
quently rich in 6-jets), it is crucial to be able to separate the c and b jets with at least 
moderate efficiency and purity. Following Ref. |75] , we examined many strategies to 
obtain this separation in the plane formed by the variables Ip-p^] and AR{fi,j) but 
without any success. The difference between the situation at the Fermilab Tevatron, 
where this strategy appears to be moderately successful, and the LHC is the kine- 
matics of the events. In contrast to the Tevatron, where jets with E^ > 25 GeV are 
readily detectable, at the LHC we have required E^ij) > 50 GeV in order not to be 
overwhelmed by mini-jet production. For this harder jet kinematics, the difference 
between mf, and rric appears to be too small to yield significant separation between 
c- and 6-jets that are not vertex-tagged. The larger contamination from 6-jets at 
the LHC only exacerbates this situation. 
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Before closing this section, we also mention one other (also unsuccessful) strategy 
that we tried for c-tagging. The idea was to utilize the difference in the distributions 
of z = E^/Ec for muons of a fixed sign of the charge from 6 or c decays. While 
the expected distributions from the quark decays are indeed significantly different, 
this strategy also fails because these quarks hadronize before they decay, and the 
^-distributions of the muons from the corresponding bottom or charm meson decays 
are essentially the same. 



3.9 Summary 

Summarizing, we have found that the use of 6-tagging enhances the SUSY reach 
of the LHC by up to 20% in a variety of well-motivated models, with the largest 
increase in reach occurring in models where rriq ^ rrig. We note that the LHC has 
used h tagging to establish regions of exclusion after 35 pb^^ of data analyzed in the 



mSUGRA/CMSSM models, as shown in Fig. [3^ from the ATLAS Collaboration 
[77j . We have also examined i-tagging, since this would have the potential to reduce 
SM backgrounds more effectively than with 6-tagging, but due to low efficiencies we 
did not obtain an increase in reach over our efforts with ^-tagging. We also attempted 
to separate the signal of third generation squarks from both the SM background and 
the SUSY signal from all other sources. Although this can be readily achieved with 
respect to the SM background, it proved to be more difficult to discriminate between 
the third generation squarks and all the other SUSY sources. 
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Figure 3.2: The statistical significance for the three NUHM model lines introduced 
in the text, assuming an integrated luminosity of 100 fb~^. All the model lines have 
^0 = and ^ > 0, with (a) > 0,tan/3 = 10, mo = mi/2, (b) > 0,tan/3 = 
10, mo = 3mi/2) and (c) m,^ < 0,tan/3 = 20, mo = 5mi/2- The solid (red) line is 
for the signal with no requirement on 6-tagging, the dashed (black) line is with the 
requirement of at least one tagged 6-jet, and the dotted (blue) line is with at least 
two tagged 6-jets. The signal is observable if the statistical significance is above the 
horizontal line at A'signai/\/-/Vback = 5. 
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Figure 3.3: The statistical significance for the LM3DM model line with mo = 
= 0, tan/3 = 10 and /i > 0, where M3(GUT) is adjusted to saturate the 
measured CDM relic density, assuming an integrated luminosity of 100 fb~^. The 
solid (red) line is for the signal with no requirement on 6-tagging, the dashed (black) 
line is with the requirement of at least one tagged b-jet, and the dotted (blue) line is 
with at least two tagged 6-jets. The signal is observable if the statistical significance 
is above the horizontal line at A'signai/V-^back = 5. 
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Figure 3.4: The statistical significance for the HM2DM model line with = 
0, tan/3 = 10, /X > and (a) uiq = ^^'^ (^) "'-o = -/Vf2(GUT). In both frames, 

M2(GUT) is adjusted to a positive value so as to saturate the measured CDM relic 
density, and an integrated luminosity of 100 fb^^ is assumed. The solid (red) line is 
for the signal with no requirement on 6-tagging, the dashed (black) line is with the 
requirement of at least one tagged 6-jet, and the dotted (blue) line is with at least 
two tagged 6-jets. The signal is observable if the statistical significance is above the 
horizontal line at -/Vgignai/V-^back = 5. 
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Figure 3.5: Plot of cross section vs squark mass for third generation squarks as solid 
(blue) curve, all generation squarks as dotted (purple) curve, all generation squarks 
and gluinos as dashed-dotted (curve) curve, with rriq = nig. 




Figure 3.6: Plot of vs mo for mSUGRA/CMSSM models, with an integrated 
luminosity of 35 pb~^ at the LHC with y/s = 7. The 6=jet channel is used to 
establish a region of exclusion in this parameter space. 
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CASE 1 


CASE 2 


CASE 3 


miQ (GeV) 


717 


854 


739 


mi/2 (GeV) 


306 


355 


361 


^0 (GeV) 


-1434 


-1708 


-1478 


tan/3 


47 


47 


47 


H (GeV) 


-372 


-428 


-477 


m-g (GeV) 


764 


879 


886 


niuj^ (GeV) 


966 


1127 


1070 


rrij^ (GeV) 


274 


316 


460 


ml (GeV) 


442 


559 


400 


(GeV) 


236 


279 


287 


£;f?'^^ (GeV) > 


150 


100 


150 


[ET{ji),ET{j2)\ (GeV) > 


100, 100 


100, 100 


200, 100 


Erihi) (GeV) > 


40 


40 


40 


meff (GeV) > 


500 


500 


600 


rij > 


5 


6 


4 


<73rd gen. (fb) 


120.2 


74.1 


80.6 


C^SUSY cont. (fb) 


1176.3 


590.6 


828.9 


<7SM (fb) 


432.6 


454.1 


580.4 


TUeS (GeV) < 


1000 


1000 


1000 


fSrd gen. (fb) 


47.2 


30.9 


20.5 


O'SUSY cont. (fb) 


109.5 


42.0 


40.0 


crSM (fb) 


141.7 


180.6 


155.1 


CSrd gen./ ^tot. bkg 


0.188 


0.14 


0.105 


-^signal/ V ^back 


29.8 


20.7 


14.7 



Table 3.6: The optimized cuts, along with cross sections for the signal from direct 
production of light third generation squarks, for Standard Model background, and 
for other SUSY contamination in the third generation signal (discussed in the text). 
The first four rows specify the input parameters for our three case studies while the 
next six rows specify and selected sparticle masses. The next several rows detail 
the choice of cuts from the set in Table [331 chosen to ameliorate the softer Standard 
Model background, along with cross sections for the third generation signal, for 
contamination to this signal from other SUSY sparticles, and for Standard Model 
background after these cuts. In the last six rows we show the cut "from above" 
discussed in the text along with our results for the various cross sections, the sig- 
nal to total background ratio (including SUSY contamination) and the statistical 
significance of the signal. 
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Variable 


Values 


nb > 


1 




300,400, ...,900 


meff(GeV) > 


800,900, ...,2000 


rij > 


3, 4,..., 8 



Table 3.7: The complete set of cuts examined for extraction of the SUSY signal 
with tagged t-jets. In addition to the basic cuts detailed in the text, we require that 
St > 0.1. 





CASE 1 


CASE 2 


CASE 3 


mi6 (GeV) 


1650 


1770 


1820 


m-g (GeV) 


1522 


1614 


1661 


rriuj^ (GeV) 


2108 


2255 


2319 


(GeV) 


714 


766 


792 


m^^ (GeV) 


744 


842 


876 


m^^ (GeV) 


533 


570 


589 


m^^ (GeV) 


279 


299 


309 


B{h tZi) 


0.64 


0.69 


0.70 


B{hi tWi) 


0.37 


0.31 


0.30 


Ef'^'' (GeV) > 


300 


500 


n/a 


rueS (GeV > 


1700 


800 


n/a 




8 


3 


n/a 


o"susY (fb) 


0.138 


0.108 


n/a 


fback (fb) 


0.0117 


0.0306 


n/a 


-^SUSY/\/^back 








top tag 


12.7 


6.14 


0.00 


16 


62.8 


52.5 


44.7 


26 


93.5 


64.0 


46.4 



Table 3.8: A comparison of the statistical significance of the LHC signal using top- 
tagging described in the text, for three different cases along the 50(10) model line 
(3.14), with other parameters as fixed by Eq. (3.13). The first few lines show the 
value of mi6 along with sample particle masses and branching fractions. The next 
three lines show the choice of cuts for the variables in Table 13.71 that maximizes 
the statistical significance of the top-tagged signal. The signal and SM background 
cross sections for these cuts are shown on the next two lines for the cut choice that 
leads to an observable signal with the greatest statistical significance. The last three 
rows compare the statistical significance of the signal using top-tagging with that 



obtained using 6-jet tagging discussed inSec. 3.4 



73 



Chapter 4 



Studying Neutralinos 
Bottom-Up at the LHC 



4.1 Introduction and Goals 



In this chapter we turn to the last of our projects introduced in Sec. 2.3 Since 
heavier neutrahnos are expected to be copiously produced via the cascade decays 
of squarks and gluinos at the LHC, we examine what we can deduce about their 
properties from a study of SUSY events. We take a bottom-up approach and focus 
on the dilepton mass (m//) distribution of lepton pairs (e or fi) produced via the 
decay Zj — t- Zfll, assuming that the two body decays of neutralinos are kinemati- 
cally forbidden. Since the leptons from neutralino decays always lead to opposite 
sign (OS), same flavor (SF) lepton pairs, i. e. e"'"e^ and /^^//^ pairs, the OS,SF 
dilepton mass distribution will play the central role in our analysis. 
We pick the mu distribution for our analysis for several reasons. Gluino and squark 
cascade decay signals in a wide variety of models exhibit readily observable rates 
of dilepton production. We can efficiently suppress the SM background to multijet, 
dilepton events with Ej^^^^ from gluinos and squarks with simple cuts on just jets in 
the event so as to minimize the loss of information contained in the mu distribution. 
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Cuts on the lepton pT or Ej^^^^ would impact the overall shape of the distribution 
making a fit to the Theoretical expectation difficult. OS, SF dileptons can also come 
from production of chargino pairs which typically lead to e~e^, fx^ fi^, e"*"//", e~ ix^ 
at equal rates. Since neutralinos always decay to OS, SF dileptons, we can use the 
readily constructible distribution N{e~ +e'^ + fJ-^ — e^/i^ — e^/i^) to statistically 
remove the "chargino contamination" to the neutralino signal. Modulo cuts, unlike 
energy or angular distributions, the rriu distribution is a Lorentz invariant, unaf- 
fected by the boost of the parent neutralino and so it is straightforward to extract 
from the data. 

The kinematic mass edge of the niu distribution has been examined in many stud- 
ies to extract m~ — m~ , e. g. m~ — m~ . This as a good starting point for an 

Zi Zf Z2 Zi 

attempt at reconstructing the masses of the neutralinos involved in the decay. We 
ask whether the mn distribution can provide us with more information than just 
this mass difference. For the case of Z2 decaying to Zi, more information would 
be possible if we could also extract -|- mg:^ which would inmediately give us a 
measurement of the neutralino masses involved in the decay. 

Other information we could extract from the mu distribution can be the relative 
sign of the mass eigenvalues. Previous studies [78j have focused on identifying this 
feature as a means of distinguishing between a higgsino-like versus a gaugino-like 
neutralino, which, we argue, is incorrect. We also ask whether it is possible to 
extract whether the Z-exchange or the slepton /^..fl-exchange contributions (shown 



in Fig. 4.1) dominate the neutralino decay amplitude. Most ambitiously we ask 
whether we can extract the parameters of the neutralino mass matrix that deter- 
mines the masses and mixing angles of the neutralinos as we detail next. Recall 
that, the neutral gauginos (A3 and Aq) and higgsinos {ipiiO and ipj^o) are not phys- 
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Figure 4.1: 3-body decay of the neutralino via virtual Z-boson exchange, and by 
virtual Il^r exchange 

ical particles with definite mass, but that these will mix to form the neutral mass 
eigenstates, the neutralinos. The neutralino Lagrangian density can be written as: 



neutralino 



2 ^ ^neutral ^ 



where 



-^neutral 

















,^ 




g 


9 I'd 


L V2 


V2 







g vu 
_g_M 


Ml 



and 



( 'I'hl \ 
V Ao J 

The entries in the mass matrix for the neutralinos, arise from the following sources 



1. The higgsino mass term fi from the superpotential. 

2. Ml and M2 from soft SUSY breaking contributions from the gaugino masses. 

3. The remaining off-diagonal terms arise from gaugino-higgsino-Higgs boson in- 
teractions, when the Higgs boson fields develop VEV's Vu = {hu) and Vd = (hd) 
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because electroweak symmetry is spontaneously broken. 

Before proceeding further we spell out the assumptions and ground rules for our 
bottom-up study. In order to keep the study as bottom-up as possible, we avoid 
specific constrained models, such as mSUGRA, where all MSSM masses and cou- 
plings are determined by a handful of parameters. 

• We assume the MSSM particle content; the shape of the ma distribution is 
then completely fixed by the neutralino parameters, 

(Ml , M2 , /U, tan(3, m^^ , m^-^ ) 

as shown in Sec. 14. 2[ 

• We will use only the shape of the mu distribution as the normalization depends 
on gluino and squark properties, decay branching fractions. 

• Assume that the two body decay channels for Z2 — )• ZiZ and Z2 — )• U are 
forbiddenj^ so our neutralino decays via the three-body mode, as shown in 
Fig. [43 

• We make a working technical assumption, for simplicity, that mj^ = 'm^^. 

We focus only upon what can be inferred from this data alone without combining it 
with information about other sparticles since that is likely to introduce other model 
dependence. 

^The first case leads to an mu distribution sharply peaked at Mz- For the second, 
mSUGRA case studies have shown that the mu distribution can be used to exclude 
off-shell slepton decays if for the test point the decay Z2 — )• II is accessible [8^. We 
make no representation as to whether or not this is possible, and if so, whether the 
conclusions of |80) extend to other models, and conservately regard the absence of 
two body neutralino decays as an assumption in this study. 
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Also we carry out our analysis for the design energy ^/s = 14 TeV of the LHC. We 
begin our phenomenological study by considering the mu distribution for the decay 
Z2 — >■ Zill where the neutralinos are gaugino-like as is typical in many models. For 
our case study we use the following parameters: 

< Mi,M2,fi,tanP,mj^,mj^ >=< 77,127,-911,10,211,211 > (4.1) 

where all mass parameters are in GeV units, and the gluino mass at 450 GeV and 
squark masses at 400 GeV, beyond the range of the Fermilab Tevatron. This is shown 



by the solid (blue) histogram in Fig. 4.2 We see that there is a sharp endpoint at 



— m~ = 50 GeV. The rates shown correspond to an integrated luminosity of 
5.5 fb^^. We introduce, 

sqn(my ) 

with the convention that r]ij = +1 for masses of equal sign , and r/jj = -1 for masses 



of opposite sign, then we can see from Figs. 4.2 that there is a difference in shape 



for the dilepton distributions corresponding to the different values of r]. 



The dotted (red) histogram in Fig. 4.2 shows the same distribution for the case 



where the sign of Mi is reversed for its value in (4.1 ) with all other parameters held 
fixed. In this case the mass eigenvalue of Zi flips its sign so that the relative sign 
r/12 goes from +1 to —1 for the dotted one. We see that the mass distribution is 
much softer in the 7712 = — 1 case. The difference in shape of the mu distribution 
due to the relative sign of the mass eigenvalues has been mentioned in the literature, 

[ZSIES]. 

The first of these references [78] attributed the differences to the gaugino-like versus 
higgsino-like nature of the neutralinos, while the second reference [7H] more appropi- 
ately to the relative sign of the mass eigenvalues. It would thus seem that we can 
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dilepton mass (GeV) 



Figure 4.2: The mu distribution from the decay Z2 — )• Zill at ^/s = 14 TeV pp 



colhder for MSSM parameters in (4.1 ) for the sohd (blue) histogram. For the dotted 
(red) histogram the sign of Mi is reversed, representing the cases for r] = +1 and 
r] = —1 respectively. 

rather easily extract the relative sign of the neutralino mass eigenvalues, but we will 
revisit this below. 

The remainder of this Chapter delves into details of information about neutralinos 
that can be extracted from the dilepton mass distribution and whether the ultimate 
goal of reconstructing the mass matrix is possible. 



4.2 The dilepton mass distribution in the MSSM 
4.2.1 Decay formula 

Let wp'' be the neutralino eigenvector components as defined in [T], and g' and g 
the hypercharge U{1) and SU{2) coupling constants. The five different coupling 
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constants in the decay process Z2 Zil~^l are: 

(1) , / (1) 

(2) , / (2) 

Let us define the following functions. 



X{x, y, z) = + + — 2xy — 2yz — 2zx 
A(mr, m-x , ) = 2m? — m~ — m~ 

t Zi2 ' t Zl Z2 

The contribution due to the Z-mediated process is, 

dVz ^ 47r^t/;Lg^(l + (-l + 4sin^M')'"V^^"''""'li'"'l^ 
dmii 12(27r)5 cos^ ("t-;^; — m|)24m~ 

((m|^ - m|^)2 - 2mf, + mfiim^-^ + m|^) + Gam^^m^^m^^,), 

where, a is the relative sign of Zi,Z2 mass eigenvalues, i.e. a = (— l)(^i+^2). 
Let us define the function, 

^rK") = -^A(m2,m|^,m|p - 4^A(m2,m|^,m|^) x 

(m~ — m^^im?- — m?) 



(2am-x m-x mfi \ 
A(mr, m^: , ) + —7 ^ 

In 



2am-^^m-^^m'ii 
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The contribution from the pure slepton-mediated decay is, 



dmi 



1 mi 

327r3 Ami 

Z2 



(4.7) 



where m^-^ , mj^ are left- and right-slepton masses respectively. Finally, define 



iz 



("I;-) = ^ (- A(m,-, m~ , + ml)^ \{ml , ml , ) 
+ ((mr — ml )(m7 — ml ) — am- m~ m?;) x 

I Zx' ^ I Z2 Zi Z2 

A(m^-, m~^ ^rn^^)+ ml + ^A(mf;, ml , ml ) 



A(m^-, m~^ , "1^2 ) + "i^ - y-^C'^i > > 
The mixed Z and slepton contribution, i.e. the cross term, is 



1 ' 
^2' 



(4. 



iz 



gmi 



dmii 327r3m~ (m?, — ml) cos^vF 
Z2 

{A~^A~wi2{l - 2sin2 ew)Fi^{miJ 

Finally, the decay width formula for the process Z2 — >• Zil^l~ is 

dV _ dVz ^ dTjj ^ dViz 
dmii dmii dmu dmn 

Clearly, the functional form of the Z and slepton contributions are completely 



(4.9) 



(4.10) 



different-looking. Note also, that each of the components in (4.10) contain the 
a term, where a is the relative sign of Zi, Z2 mass eigenvalues, i.e. a = ±1. It is this 
a term that causes the distribution to shift inward for the case a = —1 as seen in 



Fig. [42 



4.3 A Closer Look at the Test Case 



We continue our examination of our test case of eq. (4.11) and ask whether it 
is possible to infer whether the Z- or slepton-mediated amplitudes dominate the 
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decay. Toward this end, we show the slepton- mediated contribution of eq. (4.7), the 



Z-mediated contribution eq. (4.5) and the interference term eq. (4.9) in Fig. 4.3 



As expected for this gaugino-Uke case, the slepton exchange process dominates over 
the Z exchange, with a small contribution from the Z, / cross-term. In this plot, we 
have used -^^^ which shows the relative strengths of the different contributions to 
the mil distribution. In order to study whether it is possible to distinguish slepton 
and Z-mediated decays directly from the shape of the mn histogram we show the 

The lighter curve shows 



corresponding normalized contributions f in Fig. 



4.4 



Gaugino-type case contributions to total dr/dmil 



2M - mZ^ + mZ^ 
2m - mZ - mZ 



> 

C3 




20 30 40 

dilepton mass (GeV) 



Figure 4.3: Different contributions to dF/dm^; for a gaugino-like test case from the 



components of the decay formula as in (4.10). 



the shape of the Z mediated contribution while the darker curve shows that of the 
slepton mediated contribution. We can see a clear difference between the Z and 
the slepton contribution, where the latter matches exactly the shape of the total 
rail distribution. It seems that we should be able to distinguish between the two 
exchanges just from their shapes. 
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The situation is, however, not quite so simple because this conclusion presumes 



0.05 



Slepton vs Z-boson mediated neutralino 3-body decay 
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Figure 4.4: Different contributions to dT/dmn, normalized, for a gaugino-like test 
case with same MSSM inputs. 

that we know the input model parameters. The data, however, only gives us a 
reliable measurement of m~ — . In order to infer that the shape of the mu 



distribution definitely picks out the slepton-mediated contribution, we have to ensure 
that this shape cannot be reproduced by the Z-contribution with the same values 
of 2m = m~ — m~ for all values of 2M = + m~ . Toward this end, we 
examine our original test point corresponding to gaugino-like neutralinos given by 
the parameter set. 



that has an endpoint at 56 GeV. The decay is, once again, slepton dominated and 



from the slepton mediated contribution also coincides with the total. The problem, 



< Ml, M2, n, tanP, mi , m,- >=< 55, 111, -911, 10, 234, 213 > 




the corresponding 



1 dr 

r dmit 



distribution is shown in Fig. 4.5 As expected the shape 
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however, is that this shape is also the same as that obtained assuming the decay 
occurs via the Z mediated amplitude but for (m, M) = (56.38, 72.92) as shown by 
the grey line in the figure which has the same endpoint. It is clear that without any 




Figure 4.5: Normalized ran distributions showing the equivalence between Z-boson 
exchange and Il,r exchange. The values of the input MSSM parameters M\ and 
M2 are different for the respective curves, and the corresponding values of the pa- 
rameters (m, M) are shown. 

additional information, we would not be able to certify whether the shape of the 
mil histogram was the result of a Z dominated, or a slepton dominated situation. 
This observation motivated to analyze our theoretical results further. We accom- 
plish this in two ways. 

First, we look closer at our decay formula, but this time in terms of the new vari- 
ables we have introduced, with the assumption that <^ M^. Second, we look 
at our Lagrangian terms directly, for the Z and slepton contributions, applying a 
Fierz transformation, assuming that momentum transfer is much smaller than the 
virtual exchange masses. From this analysis we hope to explain the observed sim- 
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ilarity in the two contributions, observed in Fig. 4.5 The results of these analysis 
are discussed next. 



4.3.1 A Closer Look at the Decay Formula 

Consider the following transformations, in terms of new variables, m and M where. 



M + m M — m 
^Z2 = 2 '""^i" 2 ^ 



M = m2;_^-\- m^^^rn = m^;^— m^^ (4-13) 

Now apply this transformation to from the decay formula result. The question 
is, can we extract some information regarding M from the decay process. If the 
answer is affirmative, then together with the mass endpoint value, which corresponds 



to m, we would know the values for m^^ and ni^^. Rewriting our formula for 
and replacing the Z propagator with we find to leading order that 

\ = {M^-ml){rr?-ml) 

Ai=M(m^-mF,)ix[l-^]^A|[l ^ 



dV 
dm,. 



2A'P' " ^ 2A'P' 
A = 2m? - - = Ao[l - (l!^^"^] 

Ao = 2m? - ^M^ 
u I 2 

The M completely factors out in this approximation and each of the three contribu- 
tions to -j^^^ in eq. (4.10) is proportional to 



\lM\\ + a)ml + \w?] (4.14) 
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We see that the shape does not depend on M when m? <C , M% and m~ <C m? . 

Z2 ' 

We can understand this more simply using Fierz transformations as we see in the 
next section. 



4.4 Fierz Transformation Applied to Decay Formula 



The Z and slepton propagators in Fig. 4.1 can be shrunk to a point so that the 



decays occur as a four point interaction if m? <^ M~ , and the slepton is much 
heavier than Z2. We define 

Wi2 = i-if'{i)''wu 



where W12 is defined in eq. (4.3). We can write the Lagrangian for the Z mediated 
exchange as: 

C oc Wi2Zi^f,ij5fZ2Z^l{a + 675)/ (4.15) 

where 6 = 9i + 62 — 1 and a = j {Stands — cotOw) and & = | {tanOw + cotOw)- 
We define 



where A~ and B~ have been defined in eq. (4.3). Then we can write the Lagrangian 



for the left and right slepton exchange as, 

= iA'~^ 1\Z^ / + iB\^ l^nA^l + h.c. (4.16) 

The mediated amplitude comes from a four point interaction 

(r^^)„Z2a • l\ (4.17) 



We can then write 



Z2a ■ Zib = Ca r^6, li'^t/l Cyi = —-ZibT'^Z2a 



where A runs over the 16 Dirac bihnears which span the set of 4x4 complex ma- 
trices. These bihnears are classified according to their properties under Lorentz 
transformations, and consist of the following: 



{^^} = {X,75,7^757^^''"} ; //,z^ = 0,l,2,3 



(4.18) 



Substituting back into eq. (4.17) we obtain 

1 



4 2 2 



(4.19) 



from which only T = 7^ and F = ij^j^ will contribute non-zero values so that 



for the II exchange amplitude arises from, 



,1 + 75 
2 



1 - 7 5 . 



(4.20) 



• while the Ir amplitude arises from. 



Z2 ■ ht 



1 + 75^ 



(4.21) 



From eq. (4.16) the II exchange amplitude will have the form 



A.~ -1 1 A^* 

Zi ^2.1^ ^Jlii^^ Z2+ ^\ I2 -^^11^^ (4.22) 



2 



Fierz transforming the above, we obtain 

A^~ iA^-i - 1 , 



mi 2 2 



(4.23) 



but, with A^~ A~ = A~^{—if^ ^A~^{iY^ ^ and after rearranging terms (4.23) be- 



comes 



A~ A 



/ Zi-i^ X 



2 m? 

[I_(_l)ei+e2]+^^[l + (_l)g2+ei] 

Z2 



(4.24) 
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Adopting the following convention, for 6i + 62 even, choose 6 odd, and viceversa, 
then we see that when 9 is odd the neutralino current exhibits an axial vector na- 
ture, while for 9 odd it exhibits a vector nature. Then, we can write (4.24) as. 



1 IA'' ]\A^ ] 

^ 2 ii)''~'Hj'^ Pl I Zi7^ (75)^ Z2 (4.25) 

2 mf 



Summarizing, we have, 



from the Z graph. 



Zi7^(75)^Z2/V (a + 675)/ 
and from the II and Ir graphs 



From the Z graph, we can actually write l^^ {a + 675) I as, 

/ 7^^ [(a + 6)Pr + (a — 6)Pl] ^ which looks exactly like the II and Ir contributions. 
We see that ( |4.25| ) has the same form as ( 4.15). Also, the contribution from the 



right slepton would have the same structure. In our decay formula, we have used 
the massless lepton limit, so chirality is conserved, making II and Ir different final 
states which cannot interfere, explaining the absence of crossterms. 



4.5 SM background subtraction 

The real data will be contaminated by the SM background, so we still need to im- 
plement cuts to effectively remove it, while at the same time, make sure our cuts 

88 



do not alter the shape of our mn distribution significantly. Effectively this implies 
that we should focus on making cuts on jets in the event, avoiding significant cuts 
on Ej^^^^ and of course on the leptons. We use the sample gaugino-like point with 



MSSM parameters given by eq. (4.1) and show the effectiveness of our cuts on the 
SM background. We require nieptons = 2, rijets 

> 4, Ef'''' > 50 GeV and hr > 550 
GeV, where hx = ^iPTjetii = 1,4. The rriu distribution for the SUSY signal and 
the SM background is shown in Fig. 4.6P We see that the cuts effectively eliminate 
SM backgrounds. 

The question is, are there any sources of contamination to our signal still present 



Gaugino-type case w/6 GeV lepton cut & SM cuts 




Leptons ; same flavor/opp.sgn + mix flavor/opp.sgn 

susy signal 

SM b/g 

(Lum = 1 1 fb') 



60 80 100 120 140 160 180 

dilepton mass (GeV) 



Figure 4.6: Histogram showing total mu distribution of OS, SF dilepton events from 
all SUSY events (larger solid blue), together with that from SM sources (smaller 
solid red) at a \/s = 14 TeV pp collider after the cuts discussed in the text. For 
MSSM input parameters given by {Mi, M2, iJ-,tan{/3),mi,mg,mq) the gaugino-like 
case shown has MSSM parameters (77, 127, -911, 10, 211,441,441). 



at this point. Recall that our focus is with leptons of same flavour and opposite 

^We have checked that Ej^^^^ > 50 GeV cut does not significantly distort the 
shape of the mu distribution. 
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charge, as expected from neutralino decays (e~e"'" + /U"//"*") but there will be ad- 
ditional leptons from charginos in all combinations of flavor {e~e'^, fi^fi'^, e^/i^, 
e^/i"*") in equal amounts (statistically). This is seen as the long tail beyond ma = 50 
GeV, the kinematic endpoint of the distribution for the neutralino decays. In order 
to remove the contamination from charginos, we take the leptons of opposite fla- 
vor {e~^fi~+e~fj,'^) and subtract them from the total, thus, statistically eliminating 
the contribution from the charginos to the dilepton output. We show the flavour 
subtracted dilepton distribution with the corresponding SM background after cuts 
in Fig. |4.7[ We see that the tail beyond the kinematic endpoint is very efficiently 




20 40 60 80 100 120 140 160 180 

dilepton mass (GeV) 



Figure 4.7: Histogram showing the mu distribution of OS, SF dilepton events from 
all SUSY events (larger solid blue) from all SUSY events after flavour subtraction, 
together with that from SM sources (smaller solid red) also after flavour subtraction, 
at a ^/s = 14 TeV pp collider after the cuts discussed in the text.. For MSSM input 
parameters given by (Mi, M2, /U, tan{l3),m^, mg, rrig) the gaugino-like case shown has 
MSSM parameters (77,127,-911,10,211,441,441). 



eliminated. At this point, we are in a position to begin our analysis of the selected 
case studies, and we present our results in the next section. 
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4.6 Results 



To conduct our analysis we chose three representative cases, according to the com- 
position of the participating neutrahnos. 

• Gaugino-hke neutrahnos. 

• Higgsino-hke neutrahnos. 

• Mixed higgsino-gaugino-hke neutrahnos. 

We chose these cases at the beginning of 2010, prior to the LHC commencing oper- 
ations. Ah three cases have gluino and squark masses in the 400 — 450 GeV range, 
and were considered as LEP 2 and Tevatron-safe cases. The mass values for both 
gluino and squark guaranteed us a significant cross section, which for initial inte- 
grated luminosity values circa 10 fb~^ would provide us with statistically significant 
numbers of events per 3 GeV bin in our ma distributions. 

The first two case studies are presented in principle to confirm the accuracy of our 



theoretical findings presented in Sees. 4.3.1| and 4.4 where we concluded that the 



shapes of the dilepton distributions did not suffice to uniquely fit the parameters. 
This led us to consider a third case where three neutrahnos are kinematically acces- 
sible. Such a case is motivated by an interesting mechanism for obtaining CDM in 
the MSSM. The third case presented new challenges to our analysis, such as hav- 
ing three decay processes (Z3 — Z2, — Zi^andZ^ — Zi) producing dilepton 
events for the mu distribution, resulting in three different mass edges, and the LSP 
having a mixed higgsino-gaugino like nature. We focused on this case first, and as 
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it produced some interesting results with regards to our objectives, we were moti- 
vated to pursue this study further. As we were concluding this study in April 2011, 
the first results were published by the LHC, summarizing the completed analysis of 
data corresponding to an integrated luminosity of 35 pb~^. These results included 
a region of exclusion for SUSY in the (mo, mi/2) plane. Our case studies were right 
at the edge of the exclusion region at the time. Later, towards the end of August, 
beginning of September 2011, new data analysis representing integrated luminosities 
just in excess of 1 fb~^ were published, and the new regions of exclusion for SUSY 
definitely excluded our case study points. We expect that our techniques may still 
be useful for heavier gluinos and squarks in certain types of models, such as those 
having a compressed mass spectrum. Other models, such as mSUGRA, requiring 
gluinos and squarks to be at 1 Tev or more introduces a factor of 2-3 multiplying all 
the masses, so mass gaps will grow by these factors, and the separation between the 
neutralino masses will increase to a level which allows two-body decays, which will 
dominate, suppressing the three-body decays essential for our work. As we men- 
tioned before, our case studies were designed to give significant statistics at initial 
luminosities for the LHC. If we raise the gluino and squark masses beyond 1 TeV 
the cross sections will reduce significantly (by a factor of about 10), so an integrated 
luminosity of 100 fb~^, about a year of operation at design luminosity, would be 
needed in these higher mass cases. 

4.6.1 CASE 1: dominant gaugino-like neutralinos 

As mentioned above, we examine this case to verify our somewhat pessimistic con- 
clusions obtained early in this chapter, as referenced at the beginning of this section. 
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Specifically, as we try to fit our theoretical function to the data, will these fits be the 
result of the slepton mediated decay, or will it be the result of the Z-boson mediated 
decay masquerading as the slepton exchange. One way to verify this would be to 
choose the slepton mass to be very large (10 TeV) to guarantee that the slepton has 
effectively decoupled from the process, and check to see if we can obtain as good a 
fit as was obtained when the slepton was clearly involved in the process. 
Before proceeding further, a brief summary of our fitting procedure is in order. Ini- 
tially, we would like to generate an ideal, i. e. not experimentally accessible, mn 
distribution, with all dileptons being produced by neutralinos only, which we can do 
at the simulation level by identifying the parent of each lepton. In addition, we ap- 
ply only a minimal cut of Ej^^^^ = 5 GeV, no jet cuts, and no isolation requirements 
for the leptons. We require exactly two such leptons per event, of OS and SF, and 
no pt cuts on the leptons. This will give us a quasi-pure distribution, which we will 
denote by rn^f'^\ so that we can attempt to fit our theory function (;^j^)oO) to our 



expression in eq. (4.10). The shape is fit to the parameters, 

{Mi,M2,tan{/3),fi,mi) 

and an additional parameter determining the normalization. We fit by looking for 
a set of MSSM input parameters for (3^)00 and evaluating it at the midpoint of 
each 3 GeV bin, then calculating a total between the data value and the for- 
mula value for all the bin^ finally choosing the minimum of the values. The 
corresponding MSSM values plus the value for the overall constant will be our best 
fit values for the m^f^^ data. For this first case, we used rj = +1, and the MSSM 



^For the i^^ bin xf = {formula valuei — datavaluei)"^ /ni where rii is the total 
number of dileptons in bin 
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input parameter^ in eq. (4.1) where the theoretical endpoint is at 50 GeV. This 
case provides some interesting results. The number of bins is 16, and here we only 
fit one theory decay formula instead of 3 as will be required in the mixed case, so 
the number of variables is 6, (Mi, M2, /U, tan(3, 1) and the overall constant, for a net 
number of degrees of freedom = 10. Initially we attempted to identify a minimum by 
utilizing different programs available for this sort of procedure, Minuit (PAW) and 
ROOT, but it was difficult to achieve a convergent solution, and results depended 
strongly on the input of initial parameters. We proceeded by using Mathematica 
v8.1, creating a 5 parameter grid and then using a best fit. 
We include here the main results for the gaugino-like point corresponding to the fit 



for the Tn[i^^ quasi-pure case, summarized in Table 4.1 and Fig. 



4.8 



The fiat behaviour observed in the plot for higher values of m,- has a distinguishing 
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Table 4.1: results for the fit to the gaugino-like m^i'^^ case, corresponding to OS, 
SF dileptons from neutralinos w/(00) GeV lepton cuts. 



^This is the same case used to test the chargino subtraction, and SM background 
elimination. 
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Figure 4.8: \^ best fit for a gaugino-type case with leptons from neutralinos w/ 
OOGeV lepton cut, plotted vs slepton mass marginalized over other MSSM 

parameters. 

feature. In the present case, the MSSM parameters corresponding to the best-fit 
for the different values of m^- remain low, all in a relatively narrow region until 
the jump at Super- TeV slepton masses. The Ax^ remains smaller than 4 for a wide 
range of and never exceeds 8. For very heavy slepton masses we see that the Ax^ 
again attains a very low value compatible with a good fit, but for a very different set 
of gaugino masses. This is compatible with our finding earlier where we saw that it 
was not possible to distinguish between Z and slepton mediated amplitudes. 

4.6.2 CASE 2: dominant higgsino-like neutralino 

This case illustrates clearly the theoretical results obtained previously in this chap- 
ter. The mass endpoint of 24 GeV is low enough that our condition of m^^ ^ 
is satisfied over the entire range of m// values, as well as rv? <C M^. This results in 
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a complete insensitivity to the MSSM input parameters. When we create our grid 
over the MSSM parameters, we do so for a 3 GeV bin containing the value of m, 
so small perturbations from the flat minimum could occur for this reason. The 
MSSM inputs are the following: 

< Mi,M2,n,tanl3,mi^,mj^ >=< -155,170,167,10,170,170 > (4.26) 

where slepton masses are in GeV units, with m~ = 112.4 GeV and = 137.0 
GeV for a mass endpoint of 24.6 GeV. The gluino mass is 450 GeV and the squark 
masses are 400 GeV. We analyzed the m^'^'^^ case, and as expected, the behaviour 
was flat with respect to any of the MSSM parameters marginalized over the other 
parameters. Extraction of MSSM parameters is not possible. Though discouraging, 
we understand why. 

4.7 CASE 3: mixed higgsino-gaugino type neutralino 

There are SUSY models, such as the so-called High M2 DM models introduced in 
the previous chapter [79j., where there is the possibility of a Dark Matter com- 
ponent with a mixed bino-higgsino structure. In this case there will be a visible 
double mass edge from both Z2 and Z3 decays, while the third edge though present 
is not manifest. This motivated our third case study, which provided us with a 
bigger challenge for the fit, but also with very positive results. 
For this case, the gluino mass is taken at 450 GeV, while the squarks of all genera- 
tions are at 400 GeV. The remaining MSSM input parameters are, 

< Ml, M2, tanp, mi , m^- >=< -70, 400, 120, 10, 170, 170 > (4.27) 
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The resulting values for rj are, 

(^21,%l,??32) = (-1,+1,-1) 

where the slepton masses are in GeV units. In such a case, while attempting to fit 
our theory function, the fit will be sensitive to small changes in both Mi and /i , 
and to a lesser degree to changes in M2 , while exhibiting a relatively flat behaviour 
for changes in tan/3. 

The 2-body neutralino decays are kinematically suppressed, leaving dominant the 
3-body decays involving — )• Zi , Z2 — )• Zi, Z3 — )• Z2, through virtual Z , /f-slepton 
exchanges, with the relevant mass gap endpoints at: 

• m~ — m-x = 25 GeV 

^3 ^2 

• ni-x — m~ = 50 GeV 

• m-x — m-x = 75 GeV 

The corresponding dilepton mass distribution after the analysis cuts detailed ear- 
lier is shown in Fig. |4.9| for all SUSY sources together with the corresponding SM 
background. The mass edges at 50 GeV and 75 GeV are evident from the figure 
while the existence of the edge at 25 GeV can be inferred. The corresponding flavour 



subtracted distribution is shown in Fig. 4.10 The latter two endpoints are easily dis- 



tinguishable in the mu distribution, but as shown previously, we require additional 
information, and regretfully we have confirmed that + is not accessible in 
specified limiting scenarios analyzed previously. Our goal is to determine if with 
multiple endpoints and neutralinos participating, some information can be gained 
which could specify the values of the MSSM parameters of the neutralino mass ma- 
trix 
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Figure 4.9: Histogram showing total nin distribution of OS, SF dilepton events from 
all SUSY events (larger solid blue), together with that from SM sources (smaller solid 
red) at a y/s = 14 TeV pp collider after the cuts discussed in the text. For MSSM 
input parameters given by {Mi, M2, n,tan{(3),m^,mg,mg) the mixed higgsino- 
gaugino-like case shown has MSSM parameters (-70,400, 120, 10, 170,450,400). 



98 




20 40 60 80 100 120 140 160 180 

dilepton mass (GeV) 



Figure 4.10: Histogram showing the mu distribution of OS, SF dilepton events from 
all SUSY events (larger solid blue) from all SUSY events after flavour subtraction, 
together with that from SM sources (smaller solid red) also after flavour subtraction, 
at a ^/s = 14 TeV pp collider after the cuts discussed in the text.. For MSSM input 
parameters given by {Mi, M2, iJ-,tan{(3),mi,mg,mq) the mixed higgsino-gaugino- 
like case shown has MSSM parameters (-70,400,120,10,170,450,400). The solid 
(black) curve shows the result of the fit discussed in Sec. 4.7.2 of the text. 
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4.7.1 Quasi-Pure Case 



As mentioned before for the previous two case studies, we consider first the m^ii^^ dis- 
tribution for this case. We perform a least over a grid of values for (rrip Mi, M2, /i, tan/3), 
and three overall constants, each one multiplying the decay formula for the respec- 
tive neutralino decay, up to its respective endpoint. This process is more complex, 
we have a total of 25 bins, and 8 fit parameters (the SUSY paameters above plus 
independent normalizations for each of the three decays that can contribute to the 
dilepton spectrum), for a net of 17 degrees of freedom. The results from this analysis 

for the case. We have also plotted the results 



are summarized in Table 
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Table 4.2: x? results for the mixed higgsino-gaugino-like case of OS, SF leptons 
directly from neutralinos w/00 GeV lepton cuts corresponding to the fit to the 
m\^^^ distribution. 

where we have plotted vs m^ marginalized over the 
other MSSM parameters. 

The plot for the case exhibits a clear minimum with a x'^ fit of ~ 1 per degree 

of freedom (25 bins ; 8 fit variables) for a m^- ~ 200 GeV. We see the clearly rising 



for the m^'^^ case in Fig 



4.11 
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Figure 4.11: best fit to the 'm^'^^ distribution for a mixed higgsino-gaugino-like 
case, corresponding to OS, SF leptons only from neutralinos w/ (00) GeV lepton cut, 
plotted vs slepton mass marginalized over other MSSM parameters. 



on both sides of the minimum value. The best fit values are, 



{mj,Mi,M2,fi,tanl3) = (200,-82,490,126,27) 



to be compared with the input values in eq. (4.27) above. The reason that we are 
able to obtain an unambiguous fit for this case is that the ambiguity that we had 
for the single endpoint case is removed because of the contribution of the extra 
neutralinos. 



4.7.2 Quasi-Real Case 

Up to now, our analysis has been idealized in that we have pretended that exper- 
iments can identify arbitrarily soft leptons and further identify the parent of the 
lepton (since we have retained leptons from the neutralinos only). The latter issue 
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can be addressed by considering the flavour subtracted spectrum as we have abeady 
discussed. Before proceeding with results, we first discuss how we handle the issue 
of lepton cuts. 

R-function 

Our theoretical decay formula has no lepton cuts in it, while real data events will 
necessarily have cuts in the lepton px because very soft electrons and muons (< 6 
GeV in our analysis) are not readily identifiable. We expect to cause the most 
distortion to the mu distribution shape in the low mu region. Also, we expect that 
most of the effect of a not so hard lepton cut will be from kinematics rather than 
from details of the matrix element. Since we want to fit our theoretical results to the 
actual distribution we need to find a way to make an appropiate correction to our 
formula that will simulate the effect of the lepton cuts. Ideally, when calculating the 
formula, changing our integration limits for Et{1) to include these cuts would work, 
but the integrals became too complex to be able to write an analytical expression 



similar to the one we obtained for our theory function without lepton cuts, eq. (4.10 ) 
or ((ir/(im//)oo, where the subscripts indicate a GeV cut on each of the two leptons. 
In order to incorporate the effect of the lepton pT cut into our analysis, we need 



to include the effect of this cut in our theoretical fit function in eq. (4.10). Toward 
this end, we generate Monte-Carlo samples both with and without the lepton cut of 
6 GeV, and obtain the corresponding dilepton mass distributions, {dV / dmii)^^ and 
(dr/dm;i)oo- We define the ratio 

R=T^IP¥- (4-28) 

R will be larger than unity. We expect R is largest for small values of mu and 
approaches close to unity if mu is very large, for which the cut on the lepton has the 
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smallest effect. Instead of writing R as a function of mn, we write it as a function 
of the scaled variable 

mi 

X = 

— m^^ 

so that the argument of R runs between and 1. Of course the form of R will 
depend on the endpoint — m~ . If, as we expect, the effect of the lepton cut 

Z,; Z,f 

is largely kinematic R will be roughly independent of other parameters as long as 
m~ — m~ is held fixed. 

To extract the R functions for our mixed gaugino-higgsino like case with endpoints 
of 25 GeV, 50 GeV and 75 GeV, we generate four sample points having only one 
endpoint at the specified value. We can easily generate these by allowing ^ to take 
on relatively high values, so that M2 - Mi yields the desired mass gap. For the 25 
GeV endpoint we generated a gaugino like point for both the cases with Mi = it 



and also used the higgsino point with parameters given by eq. (4.26). For the 50 
GeV case We show the i?-values that we obtain from our simulation in Fig. 4.12[ 



Fig. |4.13| and Fig. |4.14| for endpoints of 25 GeV, 50 GeV and 75 GeV respectively. 
We see that the scatter for the different models is indeed small, confirming that R is 
largely determined by kinematics of the decay. The solid (red) line is our analytical 
fit for the ratios R2bi R50 and R-j^ which are parametrized as 
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„ , , 0.725* (x + 0.20) , 

„ . ^ 1.15* (x + 0.15) , 

- 1 = ,1 + (A 0.15)4 

We can then write our prediction for the spectrum of mu as 

. dr _ a . dT b dV c . dV {a 'io\ 

W^'' " R25^dma'^^^^'^ ^ R^o^dma'^-^^^ ^ W,^dm;i'^-^^^^ ^ ^ 
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Fit function for mass endpoint = 25 GeV 




dilepton mass (GeV) 



Figure 4.12: Comparison of function R25 (solid red curve) to correspond- 
ing data sets (dotted blue curves). For MSSM input parameters given by 
{Ml, M2, iJ,,tan{(3),mi,mg,mq) with all masses in GeV units, the case with 
sgn{ri) = +1 corresponds to (101,126,-711,10,175,441,441). The case with 
sgn{i]) = —1 has (101,-130,-711,10,175,441,441) and the higgsino-like case has 
(-155, 170, 167, 10, 170, 450, 400) 
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Figure 4.13: Comparison of function R^q (solid red curve) to correspond- 
ing data sets (dotted blue curves) and the high gluino mass data set 
(broken red curve w/diamonds). For MSSM input parameters given 
by {Ml, M2, fJ,,tan{f3),m^,'mg,mq) with ah masses in GeV units, the 
first four cases correspond to MSSM vals (77,127,-911,10,211,441,441), 
(77,-130,-911,10,211,441,441), to (91,141,-911,10,175,441,441), and 
(91,-144,-911,10,175,441,441), while the high gluino case corresponds to 
MSSM parameters (-71, 121, -911, 10, 315, 900, 800). 
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dilepton mass (GeV) 



Figure 4.14: Comparison of function R^^ (solid red curve) to corre- 
sponding data sets (dotted blue curves). For MSSM input parameters 
given by (Afi, M2, /x, tan(/3), m^, m^) with all masses in GeV units, the 
four cases correspond to MSSM values of (69,144,-944,10,315,414,441), 
(69,-147,-944,10,315,414,441), to (77,152,-911,10,211,441,441), and 
(77, -155, -911, 10, 211, 441, 441). 
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where a, b, c are parameters determining the normahzation, and each of the ( 



correspond to [dV / dmu)Qo in eq. (4.10). We are now ready to present our results of 
the fit to the mixed point. 



Results for the quasi-real case 



Our idealized analysis of the mixed case of Sec. 4.7.1 suggests that it may be pos- 



sible to make further progress even when realistic effects are incorporated into the 



analysis. We perform an 8-parameter fit to the flavour subtracted data in Fig. 4.10 



The values of Xmm ^^^^ each value of slepton mass are shown in Table 4.3 where 



for each value of mj the fitted values of other MSSM parameters are shown. The 





m^- GeV 


Ml 


M2 




tan/3 


61.17 


105 


-02 


210 


61 


09 


37.96 


120 


-10 


250 


61 


21 


26.00 


135 


-30 


350 


75 


32 


23.71 


150 


-43 


400 


89 


20 


23.01 


157 


-46 


420 


92 


20 


21.32 


165 


-59 


405 


109 


10 


26.24 


180 


-65 


500 


110 


23 


27.64 


200 


-85 


500 


130 


23 


30.06 


265 


-140 


600 


185 


20 


30.75 


300 


-180 


650 


225 


20 


30.41 


400 


-270 


750 


315 


20 


34.83 


550 


-405 


900 


448 


29 


39.22 


700 


-271 


600 


342 


02 


40.68 


ITeV 


-460 


800 


530 


02 



Table 4.3: results for the case of leptons from all SUSY sources w/06 GeV lepton 
cuts corresponding to the fit to the m[f^^ distribution. This case corresponds to the 



quasi-real case in the text, Sec. |4.7.2| where the cuts discussed in the text have been 
applied, together with flavour subtraction. 
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corresponding Xmin value^are shown in Fig. 4.15 where we have marginahzed over 
other MSSM parameters. We see that the slepton mass is fitted to be in the interval 




Figure 4.15: best fit to the m^f^^ distribution for a mixed higgsino-gaugino-Uke 
case, corresponding leptons from all SUSY sources w/ (06) GeV lepton cut after 
flavour subtraction and cuts discussed in text, plotted vs slepton mass marginalized 
over other MSSM parameters. Results correspond to the the quasi-real m[f^^ case 
discussed in Sec. 1122] of the text. 

150 GeV < < 175 GeV {2a). Moreover, unlike the single mass edge cases in 



Sec. 4.6 the flt is unambiguous and allows us to extract the rjij values as well as the 



relative Z and slepton mediated contributions in each case. 

4.7.3 Fitting Neutralino Mass Matrix Parameters Individually 

We look at the sensitivity for each of the other MSSM input parameters. We 
present our results for our two previous cases, the m^'i^^ and the ?7i;f^\ as our 
most ideal and most real scenarios respectively. These results are contained in the 



^In calculating the xf corresponding to the i^^ bin, we divide by rij, which corre- 
sponds to total number of dileptons in bin of OS, and both same and mixed flavour. 
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accompanying figures. Fig. 4.16 Fig. 4.17, Fig. 4.18 and Fig. 4.19 for Mi, M2, /U and 
tan{{3) respectively. The results are as expected for Afi and /u, and surprisingly 




Figure 4.16: best fit vs MSSM parameter Mi, marginalized over other MSSM 
parameters for both the m\^^^ case on the left, and the rriif^^ case on the right. 




- parameter M, - parameter 



Figure 4.17: best fit vs MSSM parameter M2, marginalized over other MSSM 
paraameters for both the 'rn^'^^ case on the left, and the m^^^ case on the right. 

sensitive for M2 but with a larger range of values. Also, as expected tanf3 exhibits 
a flat behaviour with respect to x^- the quasi-realistic fits to the data, we find 
the fitted values at the (2cr) level are 

• nil = 139 GeV to 180 GeV 

• Ml = -95 to -58 GeV 

• M2 = 360 to 475 GeV 
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Figure 4.18: best fit vs MSSM parameter /i, marginalized over other MSSM 
paraameters for both the m^'^^ case on the left, and the m^^^ case on the right. 




tanp - parameler 



Figure 4.19: best fit vs MSSM parameter tan/3, marginahzed over other MSSM 



paraameters for both the m[?^^ case on the left, and the m^"^' case on the right. 



(66) 
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= 87 to 127 GeV 



while tan{P) remains undeterred, to be compared to the input values in eq. (4.27). 
Summarizing, we set out to obtain as much information as possible from the mu 
distribution, and were surprised to learn of its relative insensitivity to m~ + 
for the case with a single mass edge. Another surprise was the ability of the Z-boson 
exchange to masquerade as a slepton exchange. This degeneracy may be resoluble 
from other data. On the other hand, the results for the favorable mixed gaugino- 
higgsino case with the double mass edge were positive, supporting our original goals 
and methodology, including our technique for simulating lepton pT cuts. 
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Chapter 5 

Conclusions and Future Outlook 



In Chapt.[3]we investigated models with an inverted squark mass hierarchy to study 
how much 6-tagging or t-tagging would increase the reach at the LHC. We also stud- 
ied ways in which 6-tagging would allow us to extract the signals for third generation 
squarks from both the SM background and the SUSY contamination from all other 
sources, which became an additional background to be eliminated. 
Now that the LHC has analyzed data corresponding to an integrated luminosity 
of just in excess of 1.1 fb~^, we find that the exclusion regions have pushed the 
gluino and squark masses above 1 TeV. This makes models with an IMH appealing, 
because with high gluino and squark masses, it is still possible to obtain light third 
generation squarks between 200 — 400 GeV. In addition, the value of S which we 
use to quantify the degree of IMH inversion, would be higher than the values we 
used for our analysis. So, for a fixed light stop (ti) mass, the value of NsusY after 
cuts would be smaller for the higher mass gluinos and squarks, increasing the value 
of our observability criteria. This makes the study of this particular area of great 
importance given the results obtained by the LHC up to date. 
In Chapt. |4]our case studies are for points which are now excluded by the LHC. As 
mentioned before, it would be worth investigating whether the fitting techniques we 
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used in our study may continue to be useful for higher values of gluino and squark 
mass as long as the neutralino mass endpoints remain at values which continue to 
suppress two=body decays, allowing three-body decays, via virtual Z or slepton ex- 
changes, to dominate. This will occur in models with a compressed spectrum, that 
have recently received some attention. We aim to pursue our study along these lines. 
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